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-REMOTE  SENSING  LABORATORY 


SUMMARY 


1 


This  docunicnt  ii'pdt  ls  t ho  rosiilts  of  a study  pi'i  forinod  (ly  tho  Roinole 
Sonsinq  Ldboro  t ory  (RSI.)  nt  t lio  University  of  Konstjs,  Lawrence,  Kansas. 

The  purpose  of  ( li  i s stuily  was  to  develop  radar  iiiraqe  simulation  techni- 
ques. This  work  is  a natuial  focus  of  over  10  years  of  prior  research 
at  RSL.  It  is  built  upon  tiu'  exiserience  developed  in  a host  of  studies. 
These  studies  include  amonij  others,  determination  of  theoretical  models 
of  radar  backseat  ter  from  realistic  terrain  and  veqetation  models, 
acquisition  of  a larrte  empirical  data  base  of  radar  backscatter  measure- 
ments, and  development  and  verification  of  the  utility  of  IDEC5  (Image 
Discrimination,  Enhancement,  Combination,  and  Sampling)  for  image  data 
ex t rac t i o n and  enhanc<*men t . The  purpose  of  this  document  is  to  summarize 
those  investigations,  to  identify  siiinificant  accomplishments,  and  to  make 
recommendations  concerning  future  research  needs  and  potential  applica- 
tions of  the  results  of  this  work. 

What  has  been  develoi'ed  is  the  capability  to  create  radar  iniaqc  simu- 
lations for  the  cieneral  case.  ^he  basic  problems  have  been  identified 
and  solved.  The  next  step  is  to  develop  specific  applications.  With 
each  application  will  come  unique,  new  problems  which  must  be  solved. 

These  applications  and  the  tiroblems  they  bring  are  the  areas  where  the 
research  efforts  must  be  concentrated  in  the  future. 

The  general  nxidel  developed  in  this  study  for  radar  image  simulation 
can  be  used  in  a variety  of  applications.  It  can  be  used  to  simulate 
the  image  products  of  diffeient  radars  and  can  accept  a variety  of  sources 
of  input  ground  truth  data.  Both  SLAR  (Side-Looking  Airborne  Radai ) 
and  PP I (Plan-Position  Indicator)  radar  images  can  be  simulated.  The 
MKxlel  can  handle  equally  well  scenes  of  terrain  having  little  or  no 
relief  and  terrain  having  significant  relief.  Potential  applications  run 
the  gamut  from  guiding  remotely-piloted  vehicles  to  photo  interpreter 
(or  pre-mission)  training. 
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A very  i inpo  r t .711 1 nii  I I's  t oiu'  iictiievecl  in  this  study  was  the  use  of 
einp  i r i r a I hacfsiattec  data  lor  1 h<‘  qrouiid  truth  tarqet  reflectance  infor- 
mation required  hy  the  model.  The  results  ot  this  study  demonstrate  the 
qreat  value  of  usiiut  empirical  hackscatter  data  in  the  simulation  model. 
This  concept  reiiMives  some  of  the  previous  limitations  placed  on  simula- 
tion and  <)reatlv  ex|)ands  potential  a()p  I i ca  I i ons  . Although  the  use  of 
empirical  data  vtas  a siqnificant  or  comp  1 i shiiien  1 , the  limited  antount  of 
these  data  av.iilahle  tesjuires  theoretical  hackscatter  models  to  extend  and 
extrapolate  the  data  lt>r  cases  not  available.  Even  with  this  limitation 
the  model  developed  has  a wide  rantie  of  applications. 

The  ohjectives  ot  the  i nves t i <ia t i on  have  been  fulfilled.  The  concept 
of  radar  imaqe  simulation  has  been  developed  and  shown  to  be  a viable 
tool  for  a multitude  of  potential  applications.  Many  areas  have  been 
identified  as  a result  of  this  study  which  require  additional  research 
in  order  that  the  promise  offered  by  radar  imaqe  simulation  be  realized. 
With  siqnificant  prociress  in  these  areas,  radar  imaqe  simulation  can  be 
exploited  to  the  fullest  extent  possible  to  aid  the  United  States  Army  in 
attaining  its  ttoals  and  mission  objectives. 


PlUFACF 


Ttu'  I'lii)  I ni'<‘ I i iKi  M'siills  siiiiii'i.i  I i /cii  in  this  i 1 ii'IK'ct  t lu'  cnKiti-’ia 

tiiui  iincl  Clin t r i hill  i ons  ol  ni.iny  i nil  i v i clua  I s . Tlio  work  was  sponsoi  ed  hy 
t tu'  Enn  i nco  r i nq  I(>(>oi|  raph  ! i 1 at'oi  a t or  i es  , Tim  lliiitod  States  Army,  Fort 
Pt' I VO  i r , \/ir(|ini<),  and  was  pnrtormrd  at  the  Remote  Sensinq  Laboratory, 

The  Univ(!rsity  ol  Kansas,  lawrinut',  Kansas . Many  pi’opte  at  these 
orqan i ^a t i ons  should  be  ieco<)ni/ed  lor  their  specific  contributions  to 
the  work  reported. 
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ABSTRACT 


I 


This  (Tocunierit  sun)i)i,i  r i ?j‘s  I ho  results  of  a radar  image  simulation 
study  nerformed  at  the  Remote  Sensing  Laboratory,  The  University  of  Kansas, 

Lawrence,  Kansas.  The  work  v;as  sponsored  by  the  Engineering  Topogra- 
phic Laboratories,  The  United  States  Arniy,  Fort  Bel  voir,  Virginia.  J 

The  goal  of  this  study  i/as  to  develop  radar  image  simulation  techni- 
ques. The  purpose  of  this  document  is  to  summarize  those  investigations, 
to  identify  significant  accomplishments,  and  to  make  recommendations 
concerning  future  research  needs  and  potential  applications  of  the  results 

of  this  wo r k 


I.n  SUMMARY 


I . I I II  t t tiiliK  I ion 

(Of  Itu'  p.ist  ? yc.if.  .1  '.tudy  li.is  bi'iMi  coiuliuttnl  ,i  t t lu>  Kcnioti'  Sons  i n<i  I .ii)- 
of.ilofy  (RSL).  Tiu‘  Univci-.ily  o(  K. ins. is,  to  'lovolop  fiiil.ir  i mukio  sinuiKitinn  .iiul 
(('.iliifc  I'xff.iclion  1 Of  tin  i ipii's  tor  it'ii.iin  iiiiapiiui  i.nldrs.  This  study  wns 

conductfd  (incK*r  conirucl  to  i (u*  Enp  i neor  i np  Topoqr.iph  i r Laboratories 
(fTL),  United  St.itos  Army.  Fort  Relvoir,  Virqinia.  Tb<>  purpose  of 
this  document  is  to  summari.'e  itiose  i nves  t i <)a  t i ons  . to  identify 
siqnitic.int  accomplishments,  and  to  make  i ecommt'ndat  ions  roncerninp 
future  research  needs  and  potential  <i()p  I i ca  t i tin  s of  t hi'  results  of 
this  wo  r k . 

We  normally  take  the  expression  "imaip'''  to  mean  some  visual 
representation  of  an  object  (or  scimk')  which  is  produced  by  lenses, 
mirrors,  or  other  optical  elements.  We  are  familiar  with  such  imapes 
when  thr'y  are  fornit'd  at  w.i  ve  U'ni|  i hs  in  the  v i s i b I (to  t i\e  human  eye) 
portion  of  tlie  e I I'c  I i om.Kiiu' 1 i c spr'ctrum  - app  rox  i m.i  1 I y iiOOO  to  7000  K. 

These  wavelcripths  are  the  reriion  of  the  e I ec  t romapne  t i c.  st>ectrum  to 
which  the  human  eye  is  sensilivr'.  It  is  possible  to  form  iniaqes  of 
objects  with  e lec  t romaqne  t i c:  w.ives  of  much  lonqer  leni|th.  In  p.ir  t i cu  1 at  . 
if  tiu'se  imaqes  are  formed  with  w.i  ve  I enq  t tis  in  the  microwave  portion  of 
the  i' I ec  t romaqne  t i c spectrum,  they  art'  called  radar  imaqes.  Althouqh 
imaqes  and  intaqinq  systr'iiis  taken  f'y  or  <i[)eratinq  in  the  microwave 
aiivl  o|itical  reiiion  of  t lu'  spi'Ctrum  have  m.iny  features  and  (iroperties 
in  common,  they  will  display  s i i|n  i f I < .in  t d i f I ei  encr's  . 

Modi'rn  in  iqinr)  r.ld.ii  s ilesiqiu'd  to  oper.lti’  i ri  complex  env  i ronnuMi  t s 
are  too  expt'nsive  to  build  just  to  sei-  il  they  will  pi'tlorm  as  desireif. 

Even  il  .i  radar  is  .ivail.ible,  I Ik*  lost  ol  opei.itinq  il  m.ikes  its  use 
un  j us  I i I i .ah  I e for  ordinary  or  I r i v i .i  I .ipp  I i r at  i ons  . Tl  is  is  c'speci.illy 
trill'  lor  airborne  or  spac.eborne  r.id.irs.  Ibest'  .ind  other  I. u tors  to  lu’  disiussed 
account  for  the  increasinq  tiemand  for  the  simulation  of  radar  imaqes. 

S iniu  la  t i op  of  radar  imaqes  means  t lia  t knowledcie  of  t fie  radar  system  and 
qroiind  scene  to  he  imaqed  are  conver  (I'd  by  suitable  techniques  into  an 


if  it  were 


inwiqe  whirl)  rt'p  rcsc'n  t s who  I the  rorlor  would  htive  "seen" 
octuolly  flown.  Rodor  s i rHU I ot i ons  are  produced,  typically,  either  by 
optical  tochn  iqiic's  (an  optical  coiuputer  is  used)  or  by  digital  techniques 
(a  high-speed  digital  computer  is  used).  flependinq  upon  the  application 
for  which  a particular  radar  i wjqe  is  to  be  simulated,  the  simulation 
effort  can  be  relatively  simple  or  it  can  be  amazingly  complex 
requiring  very  detailed  knowledge  of  the  terrain  and  vast  amounts  of 
computer  time  In  general,  the  radar  image  s i iru  I a I i on  problem  is  .not 
simple. 

Radar  imagery  possesses  several  distinct  advantages  over  optical 
photography.  The  two  prime  advantages  of  radar  are  the  facts  that 
microwave  frequencies  are  affected  less  by  meteorological  conditions, 
and,  since  radar  provides  its  own  source  of  illumination,  radar  images 
can  be  made  equally  well  during  either  day  or  night.  Cloud  cover  and 
light  precipitation  do  not  prevent  the  successful  operation  of  imaging 
radars  - hence  the  name,  all-weather  surveillance.  These  two  advantages, 
alone,  reinforce  the  value  of  imaging  radars  and,  by  implication,  they 
establish  the  potential  applications  for  imaging  radars  and,  hence, 
the  projected  value  and  applications  for  simulated  radar  imagery. 

The  capability  to  create  simulated  radar  images  has  great  potential 
value  in  many  areas.  Among  these  areas  are  such  applications  as  the 
guidance  of  retiiote  1 y-p i 1 oted  vehicles,  navigation  progress  checks, 
analysis  aid  to  the  image  interpreter,  and  training.  A flight  guidance 
computer  could  be  developed  to  use  the  simulated  r,^J^r  imagery  to  fly 
(or  guide)  the  vehicle  along  a specific  path  or  a navigator  could  use 
the  simulated  imagery  for  flight  progress  checks.  The  capability  to 
create  simulated  radar  images  can  be  seen  also  to  provide  an  active 
and  viable  analysis  and  assessment  aid  to  the  image  interpreter.  It 
would  enable  the  interpreter  to  create  for  training  purposes  (or  to 
baseline  damage  assessment  - civil  or  mi  1 i ta^'y  - In  the  event  prior 
reconna i sance  imagery  had  not  been  obtained)  virtually  any  simulated 
image  at  any  frequency  and  polarization  for  which  knowledge  of  the 
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loriMili  (i|t(iurui  tiLitli  il.Uci)  w.is  .w.i  i I nh  I o . As  vnliiahU'  ns  the  intor- 
prt'tor  is,  it  is  tiiqtily  iinlikoly  thit  cMfh  i ii  t e rp  re  I c' c is  fainiiinr 
witii  tile  itiiacie  ( roddi  Is  ol  every  radar  for  every  possible  scene. 

R.ukir  iin.Kie  s i mu  1 .1 1 i on  will  prt>vi(le  the  interprelei"  with  the  cnpnbilitv 
to  make  comparisons  outside  the  realm  of  his  experience  and  will, 
therefore,  expand  his  experituire  base  thereby  increasinq  his  efficiency 
value.  Other  potential  app 1 i ca 1 i ons  for  simulated  radar  imaqery  run 
the  qamut  from  tactical  terrain  reconna i sance  to  agricultural  crop 
health,  maturity, and  moisture  content  evaluations. 

1 . 2 Project  Summary 

Several  approaches  to  the  radar  image  simulation  problem  have  been 
investigated  at  RSL.  The  first  investigation  was  conducted  to 
analyze  the  nature  of  the  simulation  problem.  The  results  of  this 
effort  were  reported  by  Parashar*.  This  initial  effort  was  partially 
successful.  Much  was  learned  in  this  investigation  and  much  insight  was 
gained  into  the  nature  and  priority  of  the  simulation  problems  and 
parameters.  The  approach  reported  in  this  first  investigation  was 
to  model  the  test  site  as  a flat  agricultural  area.  Existing  radar 
imagery  covering  the  kinds  of  agricultural  crops  present  in  the  test 
site  were  used  to  define  the  data  base  from  which  the  simulated  image 
was  to  be  constructed.  In  an  image,  the  parameter  of  importance  is 
the  relative  variation  of  lightness  (or  darkness)  from  point  to  point. 

If  we  call  the  relative  value  of  lightness  of  a particular  image  point 
the  qreytone  of  that  point,  tht-n  the  total  range  of  variation  in  an 
image  would  be  the  greyscale.  This  first  RSL  radar  image  simulation 
work  obtained  the  greyscale  data  for  each  microwave  reflectivity 
category  (formerly  called  agricultural  crop)  to  btf  simulated  from 
existing  radar  inagery.  The  attempt  undertaken  was  to  obtain  as 
much  statistical  data  as  possible  about  homogeneous  regions  of  each 


Parashar,  S.  K.  and  A.  K.  Fung,  "Simulation  of  Radar  Imago:  Garden 
City  Test  Site,  Kansas,"  TR  23^*5.  Renrate  Sensing  Laboratory, 
The  University  of  Kansas,  May,  197^. 


coteqory  and  then  use  existing  electromagnetic  backscatler  theories 
to  interpolate  and  extend  these  limited  data.  This  was  not  a success- 
ful approach  for  the  solution  ol  the  general  simulation  problem. 

But  it  did  produce  a weal  tit  of  information  and  indicated  directions  to 
be  taken  for  future  research. 

The  second  i nves t i ga t i on  was  conducted  to  demonstrate  the  appli- 
cability of  the  radar  image  simulation  model  to  a relatively  flat 
agricultural  area  and  to  solve  soim*  of  the  problems  identified  in 

2 

the  first  study.  The  results  of  this  study  were  reported  by  Parashar  . 

The  simulated  image  was  compared  to  a real  radar  image  of  the  same 
scene.  A number  of  points  of  dissimilarity  were  identified  in  this 
comparison  but,  on  the  whole,  the  simulation  looked  like  a radar 
image.  Some  of  the  points  of  dissimilarity  were  due  first  to  the 
relative  coarseness  of  the  digital  data  base  which  represented  the 
ground  truth  information  (geometrical  and  backscatter  characteristics) 
from  which  the  simulated  image  was  produced.  The  difference  between  this 
and  the  first  study  is  the  model  used  to  extrapolate  the  empirical  back- 
scatter data.  The  area  for  which  a simulated  radar  image  was  to  be  pro- 
duced was  almost  exclusively  agricultural,  thus  it  was  decided  to  use 
Clapp's  model^  to  extend  and  interpolate  the  backscatter  data  which  were 
obtained  in  the  same  way  as  in  the  first  study.  While  for  some  applications 
Clapp's  model  may  be  a reasonable  approx i ma t i on  for  the  angular  varia- 
tion of  backscatter  data,  it  is  too  simplistic  for  the  radar  image  simu- 
lation problem  and  this  caused  a second  error  to  be  introduced.  On  the 
whole  though,  this  work  showed  that  the  simulation  of  a radar  image 
for  a medium  resolution  radar  could  be  achieved  (at  least  for  a relatively 
flat,  agricultural  area)  and  identified  more  problems  that  needed  to  be 
handled  before  the  simulation  mndelcould  be  applied  to  more  complex  areas. 
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tory, The  University  of  Kansas,  February,  1975. 
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The  next  i nvt;s  t i qa  t i on  was  conducted  to  demonstrate  the  application 

of  the  radar  itnaqe  simulation  model  to  a surface  which  has  significant 

relief.  Also,  the  microwave  reflectivity  data  used  to  calculate  the 

image  greytones  present  in  the  simulations  produced  for  the  selected 

A 5 

test  site  area  were  taken  from  applicable  scattering  theories  ’ . The 
results  of  this  study  were  reported  by  Parashar^.  The  limiting  factor 
that  prevented  more  precise  comparisons  between  simulated  and  real 
radar  imagery  of  the  same  site  was  the  specification  of  ground  truth 
information.  The  data  base  of  ground  truth  information  that  was 
prepared  for  this  test  site  was  made  by  hand  and,  consequently,  it 
was  a rough  approximation  to  that  area.  The  quality  of  the  simulation 
effort  at  this  time  was  superior  to  the  data  base.  A clear  conclusion 
to  be  drawn  from  this  study  is  that  we  need  to  define  the  minimum 
amount  of  information  to  be  included  in  the  ground  truth  data  base. 

A way  to  approach  this  probliMii  is  to  analyze  the  sensitivity  of  the 
various  parameters  which  together  make  up  the  s i mu  I a t i on  probi  em  and 
establish  for  a given  app I i ca t i on  of  radar  image  simulation  the  minimum 
required  level  of  information.  This  information  can  then  become  the 
guideline  for  building  ground  truth  data  bases.  While  on  the  subject 
it  should  be  pointed  out  that  the  most  expensive  part  of  radar  image 
simulation  is  the  building  of  ground  truth  data  bases.  The  recommenda- 
tion that  a sensitivity  analysis  be  performed  to  establish  minima 


Hevenor,  R.  A.,  "Backsca t t e r i ng  of  Electromagnetic  Waves  from  a Sur- 
face Composed  of  Two  Types  of  Surface  Roughness,"  TR  ETL-TR-71-^, 
Engineering  Topographic  Laboratories,  The  United  States  Army, 

Fort  Belvoir,  Virginia,  October,  1971. 

^ Fung,  A K.,  and  H.  L.  Chan,  "Backseat ter i ng  of  Waves  by  Composite 
Rough  Surfaces,"  IEEE  Transactions  on  Antennas  and  Propagation, 
September,  1969. 

^ Parashar,  S.  K. , A.  K.  Fung,  and  R.  K.  Moore,  "Digital  Simulation 
of  a Radar  Image  of  Pisqah  Crater  Test  Site  California," 

TR  23^-10,  Remote  Sensing  Laboratory,  The  University  of  Kansas, 
July,  1975. 
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for  qround  truth  dntn  bases  should  be  very  seriously  considered.  To 
be  able  to  simulate  a qood  quality  radar  with  a high  resolution,  a 
very  detailed  and  exact  knowledge  of  the  ground  truth  (microwave  reflec- 
tivity properties  and  geometry  of  the  site)  is  needed.  The  magnitude 
of  the  problem  is  best  illustrated  by  the  fact  that  if  it  is  desired 
to  simulate  the  image  products  of  a radar  having  a resolution  of  only 
100  feet,  the  <tround  truth  data  base  must  contain  at  least  2800 
points  per  square  mile  of  terrain  (50  feet  would  require  11,100;  25 
feet,  ^^,600;  etc).  A simulated  scene  can  easily  exceed  100  square  miles 
and  at  25  foot  reso I u t i on  t he  da t a base  would  have  to  contain  in  excess 
of  If  million  points.  The  problem  rapidly  gets  out  of  hand.  Several  other 
lessons  were  learned  in  this  investigation  that  are  worth  pointing  out  here. 
First,  more  experimental  (empirical)  microwave  reflectivity  data  needs 
to  be  collected  and  theoretical  scattering  models  need  to  be  developed 
to  support  directly  the  mission  objectives  of  the  U.S.  Army.  It  is 
strongly  beleived  that  radar  image  simulation  will  play  a very  impor- 
tant role  in  the  future  for  t!ie  U.S.  Army  (in  fact,  for  the  entire 
Department  of  Defense).  If  this  is  true,  it  will  become  increasingly 
important  to  obtain  these  data.  Secemd , it  is  necessary  to  define 
objective  image  quality  measurement  criteria.  As  the  tool  of  radar 
image  simulation  begins  to  be  used  in  more  applications  it  will  become 
increasingly  imporlant  to  be  able  to  measure  in  a quantitative  and 
objective  manner  the  quality  of  radar  image  simulations  and  their 
suitability  to  fulfill  the  task  for  which  they  were  produced. 

The  final  investigation  conducted  on  this  contract  had  the  widest 
scope  of  all  the  investigations.  This  study  was  conducted  to  investi- 
gate the  problems  laised  in  the  first  three  studies;  it  explored 
new  directions  for  the  simulation  problem  and  it  completely  generalized 
the  simulation  problem  ( t h.U  is  the  radar  image  simulation  problem  was 
viewed  as  having  iieneral,  broad-range  app I icat ions , and  the  simulation 
model  should  he  designed  to  hatulle  any  of  them  instead  of  just  one  of 
them).  The  tesulls  of  this  investigation  were  reported  by  Holtzman,  et 
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The  (H'lu'iol  nuHlel  (oi  f.idor  i ino(|c  s i rmi  I n t i cin  developed  in  (his 
study  con  be  used  to  simul.ite  the  i iiio(|e  products  of  different  radars 
and  can  accefX  a variety  of  sources  of  input  pround  truth  data.  Both 
SLAB  (S i de-Look i up  Airborne  Radar)  and  PPI  (Plan  Position  Indicator) 
radar  imapes  can  be  simulated.  The  model  can  handle  equally  well  terrain 
havinp  little  or  no  relict  ami  terrain  havin<i  sipnificant  relief.  A 
very  important  milestone  achieved  in  this  study  was  the  use  of  empiri- 
cal backscatter  data  for  the  microwave  reflectivity  data  used  in  the 
simulation.  This  study  repiesents  the  first  time  at  RSL  that  empirical 
backscatter  data  bases  were  used  in  the  simulated  radar  imape  preytone 
calculation.  Previously,  t he. backseat  ter  data  had  been  taken  by  suit- 
able statistical  technipLjes  from  a real  radar  imape  of  the  scene  to  be 
simulated.  In  this  study,  the  backscatter  data  were  taken  from  an  ex- 
tensive empirical  calibrated  data  bank  which  is  be i np  developed  at  the 
RSL  on  a contract  sponsored  by  another  apency.  The  results  of  this  study 
demonstrate  the  great  value  of  using  empirical  backscatter  data  bases  as 
the  input  microwave  data  for  the  radar  imape  simulation  problem  and,  further, 
denxjnstrates  the  potential  requirement  for  a similar  effort  to  be  conducted  on 
behalf  of  the  qoals  and  mission  objectives  of  the  Army.  No  new  data 
bases  were  constructed  in  this  study  and  the  simulation  model  was 
vastly  superior  to  the  available  pround  truth  data.  A very  pood  data 
base  needs  to  be  developed  so  that  the  true  fidelity  of  the  simulation 
model  can  be  appraised.  All  of  the  first-order  effects  (e.p.  the 
niost  important  effects)  were  incorporated  in  the  simulation  model. 

These  include  shadow,  layover,  near-range  compression,  fading,  etc. 

Not  all  of  the  effects  were  shown  in  the  sample  radar  image  simula- 
tions presented  in  this  last  report  bacause  the  data  base  was  unsuit- 
able; it  was  flat  terrain.  Rut  the  problems  were  solved  and  applications 
will  be  shown  in  this  report.  An  alternative  to  the  static  simulation 
model  was  explored  in  which  a human  being  (photo  analyst,  etc.)  was 
placed  in  the  decision  process  of  the  computer  software  in  an  interactive. 


Holtzman.  J.  C.,  V.  H.  Kaupp,  R.  L.  Martin,  C.  E.  Komp,  and  V.  S. 

Frost,  "Radar  Imape  Simulation  Project:  Development  of  a General 
Simulation  Model  and  an  Interactive  Simulation  Model,  and  Sample 
Results,"  TR  23^"13.  Remote  Sensing  Laboratory,  The  University 
of  Kansas,  February,  1976. 
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if.)l-tinio  .1'.  tin-  proi|toin  r.in)  ope  ra  t i on.i  I mode.  This  interactive 

simulation  teclinii|ue  was  slu>wn  to  he  hiuhiv  flexible  and  t«)  be  a siiini- 
I leant  i mt)i  ovemen  I ovi'i  i onven  I i oua  I t ih  Itn  i cpies  which  reiiuiie  .uIIum - 
once  to  inflexible  formats  with  luird-codino  of  the  sii’iulation  instructions 
prior  to  production,  and  thus  the  simulation  parameters  and  scene  cate- 
^^ori^ations  are  difficult  or  im()ossihle  to  alter  easily  and  rapidly. 

It  is  felt  that  t tie  value  oi  interaction  between  computer  and  analyst 
has  been  stiown  in  this  s turfy.  By  i nc  1 ud  i net  the  human  in  the  decision 
process. a wfiole  greater  tlian  the  sum  of  its  parts  is  created. 

The  computer  is  a vr'ty  useful  tool  hut  its  dec  i s i on-mak  i no  capabilities 
are  very  limited.  By  usitK]  t tie  computer  to  maximize  its  strengths 
(e.q.  rapid  manipulation  of  vast  anxiunts  of  data)  and  using  the  trained 
human  to  make  dr'cisions,  tremenrlous  improvements  can  be  achieved  in 
the  radar  image  simulation  piohlem.  By  extrapolating  this  capability 
to  other  areas  it  seems  natural  to  infer  that  the  interactive  mode 
is  the  obvious  rirethod  of  choice  for  preparation  of  data  bases.  Much 
of  the  problem  of  riafa  base  coir^  t rue  t i on  is  centered  around  the  drudge 
of  data  manipulation;  the  computer  is  very  good  at  this  task.  On  the 
other  hand,  many  decisions  must  be  made  regarding  the  nature  of  the 
terrain  from  which  a data  base  is  to  be  constructed;  humans  are  very 
(lood  at  this  task  but  computers  arc  not.  The  optimum  method  for 
data  base  construction  would  combine  the  strengths  of  both  computer 
and  human;  t lie  computer  for  data  manipulation,  image  enhancement 
techniques,  statistical  image  classification  and  pattern  recognition 
fechnigues,  .ind  lire  lium.rn  for  his  own  special  expertise:  applyirui  his 
yt'ars  of  Itainirui  and  praclirc  to  make  imi>ortant  derisions  regard!  nr) 
t fie  intrinsic  p i <1)10 1 1 i es  o I I lie  ter  r a i n 

What  has  been  rU'veloiied  is  the  capability  to  err'ate  radar  image 
simulations  for  a host  ol  app  I i 1 a r i ons  . Ttie  basic  iirohlerns  have  been 
identific'd  and  solved.  Wh.it  rem.iins  is  to  develop  the  app  I i ca  t i ons  . 

With  each  atiplication  will  comt'  unir]ue,  new.  problems  which  must  be 
solved.  This  is  tin-  area  whc're  much  work  remains  to  be  done. 


H 


I . 3 Potential  Applications  of  Radar  Image  Siniulation 


The  radar  imaqe  simulation  model  developed  on  this  contract  is 
a rigorous  mathematical  model  which  treats  the  radar,  scene,  and  imaqe 
medium  as  a closed  system.  The  variables  which  comprise  the  model 
of  this  closed  systeni  can  be  specified  as  precisely  as  is  desired.  The 
application  for  which  the  simulated  radar  imaqe  is  to  be  created  will 
dictate  the  degree  of  detail  that  is  required  for  each  of  the  simula- 
i tion  parameters  and  the  data  base. 

: Since  this  simulation  model  is  a rigorous  mathematical  model, 

its  potential  applications  are  practically  limitless.  Rather  than 
list  and  discuss  the  applications  which  immediately  come  to  mind  when 
thinking  about  radar  image  simulation  (e.g.,  navigation,  terrain  analysis, 
\ damage  assessment,  guidance,  image  interpretation,  feature  extraction, 

mission  planning,  terrain  classification,  training,  etc.)  we  wish  to 
discuss  a few  less  obvious  potential  app 1 i ca t i ons . For  instance,  radar 
image  simulation  could  become  an  invaluable  tool  in  future,  radar  system 
development  programs.  S i mu  1 a t i on  cou  1 d be  used  to  determine  the  optimum 
frequency,  polarisation,  antenna,  transmitter  characteristics,  receiver 
characteristics,  and  etc.,  for  discrimination  of  the  desired  targets. 

In  other  words,  simulation  could  be  used  to  determine  the  minumum 
specifications  and  tolerances  for  the  various  (or  just  the  critical) 
sensor  parameters.  It  could  also  be  used  to  evaluate  the  testing 
criteria  themselves  to  determine  the  ability  of  the  tests  to 
measure  the  desired  effects. 

Once  candidate  radar  systems  are  designed,  simulation  could  be 
used  to  evaluate  and  compare  competing  designs.  Their  applicability 
to  specific  program  goals  and  their  capabilities  to  meet  those 
objectives  could  be  evaluated.  Simulation  could  be  used  to  compare 
the  capabilities  of  existing  radars  and  measure  their  applicability 
to  the  goals  and  objectives  of  a particular  program  and.  thus,  could  be 
used  to  help  make  the  decisions  of  whether  to  build  a new  system  or 
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to  use  an  exist  io(|  sysleni.  Simulation  could  be  used  to  evaluate  the 
parameter  tradeotts  and  to  help  determine  the  radar  specifications  for 
a new  program.  SiiMulation  could  also  be  used  to  define  and  evaluate 
scientific  as  well  as  eng i tieer i ng  testing  criteria.  It  could  be  used 
to  produce  sample  data  for  evaluating  the  testing  criteria,  data 
handling  systems,  and  for  training  personnel. 

In  other  words,  precisely  because  this  model  is  mathematically 
rigorous  it  is  extremely  flexible.  This  flexibility  and  utility 
should  not  be  lost  but  should,  instead,  be  exploited  to  the  ful lest 
extent  possible  to  help  the  U.S.  Army  attain  its  goals  and  mission 
object i ves . 
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i . 0 KAD AK  SIMULATION  MODI  L 
?.l  Simulation  Introduction 

An  an.i  I y t i i, .)  I aiipi  i cu  li  1 1 ■ I lio  p i oli  I oni  o I t I lo  •.  i mu  I a I i on  o I i ada i 
imaqery  was  developed.  This  approach  resulted  in  a unique  description  of 
the  radar  image  simulation  problem  by  considering  the  image  product, 
radar  system  and  ground  scene  to  be  a closed  system.  Figure  1 is  a 
conceptual  model  for  this  closed  system.  Reference  to  this  figure  will 
show  the  relationships  between  the  three  basic  parts  of  the  closed 
system:  image  medium,  radar  system,  and  ground  scene.  A simplistic 
explanation  of  this  system  would  be  as  follows.  A transmitter  generates 
short  bursts  (or  pulses)  of  energy  (at  microwave  frequencies)  which  are 
blocks  of  energy  traveling  at  the  velocity  of  light  propagated  into 
space  by  means  of  a directional  antenna.  The  energy  is  confined  to 
illuminate  the  surface  of  the  earth  in  a narrow  path.  At  any  one 
instant,  the  area  of  the  earth's  surface  illuminated  by  the  transmitted 
pulse  is  limited  in  the  direction  of  the  pulse's  propagation  (range 
direction)  by  the  physical  length  of  the  transmitted  pulse  and  in  the  ortho- 
gonal direction  (azimuth  direction)  by  the  beamwid th, rea I or  synthetic, of 
the  directional  antenna.  Objects  which  are  reflective  to  the  energy  at  a 
particular  frequency  or  wavelength  will  re-radiate  a fraction  of  the  trans- 
mitted energy  back  to  the  antenna  (this  is  called  backscatter  and  the  per- 
centage of  energy  re-radiated  back  in  a specific  direction  per  unit  area  by  a 
particular  object  is  called  its  scattering  coefficient  - o”).  This  energy 
will  be  received  by  the  antenna  and  will  be  converted  by  the  radar  receiver  to 
a video  signal.  The  magnitude  of  this  video  signal  at  any  instant  in  time 
is  determined  by  the  scattering  properties  of  the  objects  having  the  same 
round-trip  delay  time  from  transmitter  to  ground  to  receiver.  That  is, 
for  each  pulse  of  energy,  the  receiver  processes  the  received  energy 
according  to  the  time  elapsed  from  the  time  the  energy  was  transnntted; 
the  energy  re-radiated  from  objects  which  are  closer  to  the  radar  will 
be  received  and  processed  before  the  energy  from  objects  which  are  farther 
away.  This  video  signal  can  he  processed  in  a number  of  ways.  The  object 
is  to  convert  the  signal  into  an  image  and  many  ways  are  available. 


Closed  System  Model  for  Simulation  of  Imaging  Radar 


Figure  1.  Closed  System  Model  for  Simulation  of  Imaging  Radar 


1(11  example,  i(  c <in  he  recoi 'led  directly  on  sicpial  f>lm;  it  can  he 
converted  to  dig  it.)  I data  and  stmeil  on  marinelic  tape;  it  can  he  displayed 
by  intensity  modulatinq  the  electron  beam  of  a CRT  (Cathode  Ray  Tube), 
and  plnrl  o(|  I .rphcnl  , or  it  can  he  ptixessi'd  by  othiM,  mori'  elaborate  methods. 
The  I i na  I i m<)(ie  (the  v i sn.i  I record  o(  the  radar  return  si  (trial)  depends  on 
variations  in  the  relative  strengths  of  the  sirinal  returned  from  the 
dillerenl  |iatt'.  ol  the  are.r  imaiied  to  prruluce  contrasts,  erhies  , and  the 
range  of  brightness  in  the  image.  The  tvxj  primary  factors  deterrni  n i ng  the 
strength  of  the  sirinal  ribserved  on  the  final  image  produced,  and  conse- 
quently the  brightness  of  an  image  point,  are  qeorretric  and  dielectric  pro- 
perties of  the  target. 

2 . 2 Simulation  Mode  1 

The  development  of  the  radar  image  simulation  model  will  not  be 
repeated  here.  The  development  is  reported  by  Holtzman,  et  al.^  and  the 
results  are  summarized.  Recalling  that  the  relative  brightness  between 
white  and  black  of  a point  in  an  image  is  called  the  greytone  of  that 
point,  then  the  parameter  we  wish  to  calculate  for  each  point  in  an  image 
is  the  greytone.  The  greytone  for  each  point  in  an  image  can  be  computed 
from  the  following  equation; 


Gp  = The  instantaneous  greytone,  including  fading,  to  be  calculatad 
^2  for  each  discrete  point  for  scatterers  belonging  to  backscatter 

category  ' c ' ; 

G = The  greytone  value  added  to  the  value  computed  for  each  point 

1 to  calihi.ilo  I he  i .ini|i'  ol  hi  i (|hl  in  an  im.np-  .\(ioidiiii| 

to  a known  reference; 

Y = A property  of  the  image  medium  (in  this  case,  film,  i.e.  image 
iiH'd  i urn  transfer  function); 

P.p  = The  transmitter  output  power  of  the  radar  to  be  simulated; 

2 


^Holtzman,  J.  C.,  V.  H.  Kaupp,  R.  L.  Martin,  E.  E.  Komp , and  V.  S.  Frost, 
"Radar  Image  Simulation  Project:  Development  of  a General  Simulation 
Model  and  an  Interactive  Simulation  Model,  and  Sample  Results,"  TR 
23^*13,  Remote  Sensing  Laboratory,  The  University  of  Kansas,  Feb.,  1976 
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^2  ~ the  iiround  sfiot  resolution  cell  iMuniinatea  Dy 

the  radar  to  be  simulated; 

AA|  = Area  of  the  ground  spot  resolution  cell  illuminated  by  the 
cal i bra  tor ; 

G2  = One-way  gain  of  the  antenna  of  the  radar  to  be  simulated  (in 
direct! on  of  ' A^ ) ; 

G|  = One-way  gain  of  the  antenna  of  the  calibrator  (in  direction 


of  \Aj); 

^2  = Wavelength  of  the  electromagnetic  energy  transmitted  by 
the  radar  to  be  simulated; 

I,  j = Wavelength  of  the  electromagnetic  energy  transmitted  by  the 
ca  I ; bra  tor ; 

R2  = The  distance  from  the  antenna  of  the  radar  to  be  simulated  to  each 
ground  resolution  cell; 

Rj  = The  distance  from  the  antenna  of  the  calibration  system  to  the 
reference  ground  spot; 

k2.kj  = Constants  which  depend  upon  the  exposure  time  and  on  the 
film  processing  and  development; 

~ Constants  of  proportionality  relating  the  return  electromagnetic 
power  received  by  the  antenna  to  intensity  on  the  film; 

N = The  number  of  independent  samples  contained  in  uncorrelated 
resolut ion  cells; 

RN  = A Rayleigh  distributed  random  number  between  zero  and  one 
having  zero  mean  value  and  variance  equal  to  one. 

Not  explicitly  shown  in  this  equation  (I)  are  the  methods  by  which 
the  various  properties  of  the  imaging  problem  have  been  satisfied.  For 
instance,  since  the  earth  is  not  smooth  and  flat,  when  calculating  the 
angle  between  a line  perpendicular  to  the  surface  of  the  earth  at  a point 
and  a line  from  that  point  to  the  antenna  of  the  radar,  the  components  of 
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that  the  model  will  ho  invalid  in  another  configuration.  The  model  is 
completely  general  and  rigorous.  The  model  should  be  simplified  accord- 
ing to  the  specific  regui foments  of  each  application.  Only  by  doing  this 
will  the  true  value  of  the  model  he  realized.  Unfortunately,  it  is  not 
known  towhat  level  of  detail  it  is  necessary  to  specify  the  ground  truth 
data  for  any  particular  application.  Since  this  represents  the  major 
cost  in  radar  imago  simulation  (()round  truth  data  bases)  the  problem  needs 
to  be  studied  and  the  minii.um  degree-of -spec i f i ca t i on  of  the  ground  truth 
data  bases  for  specific  applications  needs  tf)  be  established.  Only  after 
this  is  accomplished  can  the  true  simplicity  and  value  of  the  radar  image 
simulation  model  developed  here  be  realized. 

This  model  has  been  tested  atjainst  a data  base  consisting  of  geomet- 
ric solids.  The  data  base  of  geometric  solids  is  a test  data  base  that 

was  constructed  to  evaluate  the  output  responses  of  the  simulation  model. 

9 

The  geometric  data  base  has  been  reported  and  is  included  as  Appendix  B 
to  this  report.  The  geometric  data  base  has  been  an  invaluable  testing 
tool  because  it  provides  us  with  the  capability  to  measure  objectively  the 
response  produced  by  our  simulation  model;  the  input  data  are  completely 
known  and  the  output  response  is  completely  predictable.  Precisely 
because  the  output  response  is  calculable  can  this  test  tool  be  used  to 
measure  the  validity  of  our  simulation  nx3de  I . Sample  results  of  the 
application  of  our  simulation  model  to  the  geometric  data  base  are 
included  in  Appendix  A. 

2.3  Si  mu  1 a t ion  Cone  I us i ons 

At  this  stage  in  the  development,  the  radar  image  simulation  model 
is  superior  to  the  available  data  bases  (at  RSL)  of  real  scenes  for  which 
we  can  test  it.  The  next  step  will  be  to  compare  a simulated  radar  image 
and  a real  radar  image  of  the  same  site.  A very  good  data  base  is  in  the 
process  of  being  built  to  serve  this  purpose.  Upon  completion  of  that 
data  base,  the  next  stage  of  refinement  regui red  to  make  this  simulation 
model  applicable  to  the  various  potential  uses  for  simulated  radar 
imagery  will  be  determined.  It  should  be  said  that  by  further  refinement 


9 

Komp , E.  0.,  V.  H.  Kaupp,  and  J.  C.  Holtzman,  "Construction  of  a Geometric 
Data  Base  for  Radar  Image  Simulation  Studies,"  TR  319-1.  Remote  Sens- 
ing Laboratory,  The  University  of  Kansas,  July,  1976- 

16 


is  not  necessarily  mean'  adilitiotial  complexity.  In  fact,  simplification 
of  the  simulation  model  dat<)  requirements  is  souqht.  Of  the  three 
components  of  the  clos('d  system  (imaije  product,  rada  system,  and  qround 
scene)  which  are  included  in  this  model,  the  area  In  which  simplification 
would  have  the  greatest  impact  is  in  the  ground  scene  data  base.  Ground 
scene  data  bases  are  terribly  expensive  to  make  and  any  reduction  which 
can  be  made  in  the  minimum  level  of  detail  required  in  the  specifica- 
tion of  the  terrain  data  translates  into  dollars  saved.  A study  needs  to 
be  conducted  to  establish  these  criteria  for  each  of  the  potential 
applications  for  radar  image  simulations. 


} . 0 CO NCLUSIONj  an ni_ R E£0 mm cndat  i ons 
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Whnt  lias  been  lievelopeH  is  the  capability  to  create  radar  image  simula- 
tions fcir  a wide  variety  of  applications.  The  basic  geometrical  problems 
have  been  identified  and  solved.  What  remains  is  to  develop  the  applica- 
tions. With  each  application  will  come  unigue,  new  problems  which  must  be 
solved.  These  applications  and  tlie  problems  they  bring  are  the  areas  where 
research  effort  must  be  concentrated  in  the  future. 

The  study  was  conducted  to  investigate  the  applicability  of  radar 
image  simulation  to  assist  the  U.S.  Army  in  attaining  several  specific  mis- 
sion objectives.  Each  phase  of  tlie  study  explored  new  directions  and 
broadened  the  general  applicability  of  the  concept  (radar  image  simulation) 
to  more  areas.  As  the  investigations  continued,  it  became  apparent  that 
radar  image  simulation  would  in  the  future  become  an  increasingly  important 
tool.  Radar  image  simulation  has  not  lived  up  to  its  potential  in  the  past 
because  of  several  very  real  limitations.  These  limitations  are  being 
overcome  with  the  use  of  high-speed  digital  computers,  with  the  develop- 
ment of  empirical  microwave  backscatter  data  banks,  with  more  automated 
ways  to  develop  the  ground  truth  data  bases.  By  the  conclusion  of  this 
study  the  basic  theoretical  model  had  been  developed.  The  basic 
technical  problems  had  been  defined.  A rigorous  mathematical  nxjdel  had  been 
constructed.  Remaining  were  the  difficult  problems  associated  with  apply- 
ing the  model  to  particular  radar  systems  for  specific  uses.  Many  radar 
phenomena  are  implicitly  included  in  the  model  developed  but  depend  upon 
imp  1 emen ta t i on  of  the  model  for  explicit  representation.  These  phenomena 
together  with  the  problems  associated  with  constructing  ground  truth  data 
bases  almost  invariably  depend  upon  the  parameters  of  the  radar  system  being 
modeled  and  the  intended  use  of  the  resultant  simulated  imagery.  Even 
though  a rigorous  mathematical  model  has  been  developed  which  is  practical 
to  use,  much  future  research  effort  is  still  reguired.  To  maximize  the 
value  of  this  simulated  capability  this  research  effort  must  be  concen- 
trated on  the  problems  associated  with  specific  applications  of  the  simula- 
tion, ground  truth  data  base  construction,  empirical  data  acquisition 
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ond  theoretical  modlinq,  large  ii'iage  haridlitu)  problems,  and  cultural 
scene  s iniu  la  t ion . These  are  the  |)rime  areas  where  immediate  research 
is  needed. 


The  general  model  de!velo)ied  in  this  study  for  radar  image  simula- 
tion can  be  used  in  a variety  of  applications.  It  can  be  used  to 
simulate  the  image  products  of  different  radars  and  can  accept  a 
variety  of  sources  of  input  groLind  truth  data.  Both  SLAB  (Side-Looking 
Airborne  Radar),  real  and  medium  resolution  synthetic  aperture,  and  PPI 
(P I an-Pos i t i on  Indicator)  radar  images  can  be  simulated.  The  model 
can  handle  egual ly  well  scenes  of  terrain  having  little  or  no  relief 
and  terrain  having  significant  relief.  The  model  is  developed  by 
Holtzman,  et  al.^  Potential  applications  run  the  gamut  from  guiding 
remotely-piloted  vehicles  to  photo  interpreter  (or  pre-mission)  training. 

A very  important  milestone  achieved  in  this  study  was  the  use  of 
empirical  backscatter  data  for  the  ground  truth  target  reflectance 
information.  The  results  of  this  study  demonstrate  the  great  value  of 
using  empirical  backscatter  data  in  the  simulation  model.  This  concept 
removes  some  of  the  previous  limitations  placed  on  simulation. 

Although  the  use  of  empirical  data  was  a significant  accomplishment, 
the  limited  amount  of  these  data  available  reguires  theoretical  back- 
scatter models  to  extend  and  extrapolate  the  data  for  cases  not  available. 

Even  with  this  limitation  the  model  developed  has  a wide  range  of  appli- 
ca  t i ons . 

Most  terrain  can  beniodcled  as  collections  of  a number  of  homogeneous 
regions.  These  homogeneous  regions  are  typically  much  larger  than  the  resolu- 
tion element  of  the  radar  system  being  iiKideled  so  they  are  normally  called 
distributed  targets.  The  isode I for  radar  image  simulation  developed  here 
handles  distributed  targets  exceptionally  well.  Cultural  targets  (hard 
targets)  are  not  handled  so  well  by  this  method  because  of  the  extreme 
complexity  of  adeguately  specifying  the  geometric  and  dielectric  properties 
of  this  class  of  target.  Another  model  involving  optical  techniques 

1 
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reported  by  Frost  , et  al..  is  the  choice  for  simulation  of  cultural 
targets.  The  optimum  method  for  simulating  a scene  consisting  of  both 
distributed  and  cultural  targets  would  combine  the  strengths  of  both 
techniques;  digital  simulation  techniques  for  distributed  targets  and  optical 
techniques  for  cultural  targets.  Given  this  method  for  the  simulation 
model,  the  capability  to  produce  high-quality  simulated  radar  images  far 
exceeds  the  capabilities  to  produce  high-quality  data  bases. 

3 . 2 Recomtnenda  t i ons 

The  objectives  of  the  investigation  have  been  fulfilled.  The  concept 
of  radar  image  simulation  has  been  developed  and  shown  to  be  a viable 
tool  for  a multitude  of  potential  app I i ca t ions . Many  areas  have  been 
identified  as  a result  of  this  study  which  require  additional  research  in 
order  that  the  promise  offered  by  radar  image  simulation  be  realized.  With 
signi^^icant  progress  in  these  areas,  radar  image  simulation  can  be  exploited 
to  the  fullest  extent  possible  to  aid  the  U.S,  Army  in  attaining  its 
goals  and  mission  objectives.  Without  significant  progress,  radar 
image  simulation  will  continue  to  be  an  expensive  tool  with  only  limited 
app I i ca t ions . The  following  discussion  treats  3 partial  list  of  the 
most  obvious  areas  which  need  work.  The  discussions  are  based  upon  the 
problems  and  achievements  encountered  during  the  performance  of  this 
study.  Importance  has  been  assigned  to  these  problem  areas  according  to 
their  s i gn  i f i caiice  as  they  appeared  to  the  investigators. 

3.2.1  Feature  Extraction 

The  utility  and  versatility  of  radar  image  simulation  is  dependent 
upon  automating  the  problems  of  ground  truth  data  base  definition. 

Clearly,  the  biggest  single  problem  to  be  handled  is  that  of  identifying 
the  geometric  and  e I ec t romagne t i c properties  of  the  scene  which  are  to  be 
transferred  into  the  data  base.  The  term  most  often  applied  to  this  task 
is  feature  extraction. 


Frost,  V.  S,  J.  L.  Abbott,  V.  H.  Kaupp,  and  J.  C.  Holtzman,  "A  Mathe- 
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Classical  te'chniciues  for  fcaiort'  (Extraction  are  primarily  manual. 
Typically,  a photo  interpretc'r  scans  t fu*  i n te  1 1 i <ience  data  and  draws 
upon  his  interpretation  experience  to  decide  what  information  to  transfer 
manually  to  the  data  base  under  construct  ion.  These  decisions  are  made 
with  as  few  diuital  computer  imarie  enhancement  techniques  as  possible. 

This  reticence  to  use  available  enhancement  routines  is  caused,  in  part, 
by  the  very  nature  of  the  auloi'iatic  routines.  They  are  not  generally 
applicable  to  any  but  specific,  we  1 1 -st ructured , test  cases.  In  addition, 
use  of  these  techniques  requires  that  the  interpreter  also  be  a computer 
expert.  Moreover,  the  interpreter  loses  control  and  visibility  of  what 
he  is  trying  to  accomplish  v;hen  he  enters  the  computer  world  of  automatic 
land-use  classification,  or  pattern  recognition,  or  region  definition,  or 
ad  infinitum.  These  reasons  have  serious  ramifications  for  feature 
extraction  and,  consequently,  data  base  construction;  they  cost  money. 

They  cost  money  in  the  sense  that  it  takes  a much  longer  time  to  extract 
the  features  for  a data  base  than  might  otherwise  be  necessary;  data 
are  manipulated  by  hand  and  the  best  information  may  not  have  been 
obta i ned . 

Clearly,  a tremendous  improvement  of  the  product  developed,  resources 
expended,  and  time  required  could  be  obtained  if  a workable  marriage 
between  computer  and  interpreter  could  be  arranged.  The  computer  is  very 
good  at  manipulating  vast  amounts  of  data  in  short  periods  of  time;  the 
human  is  not.  The  human  is  beyond  comparision  when  it  comes  to  drawing 
upon  learning  experience  to  make  decisions.  The  computer  excels  at 
clearly  defined  repetitive  tasks,  at  statistical  analyses,  at  image 
enhancements.  A cooperative  approach  in  which  the  human  is  used  to  make 
decisions  and  guide  the  processing  directionof  the  software,  and  the 
computer  is  used  to  manipulate  the  data  rapidly  and  easily  and  to  remove 
the  drudge  from  the  human  would  be  optimal  in  the  sense  of  maximiz- 
ing the  return  for  resources  expended  and  minimizing  the  time  and  effort. 
This  cooperative  approach  is  called  interactive  feature  extraction,  or 
automated  feature  extraction  (not  automatic  feature  extraction  since 
this  is  impossible  with  the  slate  of  the  art  available  today). 

The  concept  of  interactive  feature  extraction  uses  the  human  for  his 
I specific  strengths  and  the  computer  for  its  specific  strengths.  In 

interactive  feature  extraction,  the  computer  is  used  to  display,  enhance. 
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ni.in  i pii  1 .1 1 1* , .111(1  ottuTwiso  .i  i d t In*  liinn.m  interpreter  as  he  performs  his 
function.  Viewed  .mother  way,  the  tuiman  is  used  to  make  decisions  and  to 
()uide  the  cotnputer  in  real-time  as  the  programs  run.  Interaction  can  be 
accomplished  by  (liviiu)  the  interpreter  a f eiv  basic  tools  with  which  to 
communicate  his  decisions  to  the  computer;  a kt'yhoard  for  comm.ands  and 
joystick  for  direct  spec  i f ic.it  i on  .ire  probably  t fu'  minimum  to  be  provided. 
Given  these  capabilities,  the  data  h.ise  can  be  built  directly  as  the  feature 
information  is  prcicessed  and  rlecisions  are  made.  Boundaries  separatinc) 
different  regions  c.in  be  specified  directly  by  the  interprc'ter  .ind.  while 
the  human  is  analysing  the  next  problem  area,  the  coniputer  can  build  the 
symbolic  data  base  immediately  and  display  the  results.  Depending  upon  the 
level  of  sophistication  of  the  interactive  software  and  the  computer  and 
display  complex,  tremendous  savings  of  resources  and  improvements  in 
efficiency  and  guality  of  the  finished  product  are  visualized.  Given  an 
interactive  feature  extraction  system,  special  emphasis  could  be  built  in 
to  maximize  the  use  of  the  intelligence  data  normally  available  from  which 
to  define  the  geometry,  dielectric  properties,  and  elevation  data  which  are 
reguired  by  radar  image  simulation. 

3.2.2  Sensitivity  Analysis 

The  utility  and  versatility  of  radar  image  simulation  can  be  improved 
and  the  cost  reduced  if  the  minimum  level  of  detail  required  to  be  in 
the  data  base  for  specific  applications  of  radar  image  simulation  can  be 
determined.  As  previously  noted,  the  most  expensive  part  of  the  radar 
image  simulation  process  is  the  building  of  the  digital  data  base.  If 
it  can  be  determined,  for  a specific  application,  that  the  level  of  detail 
in  the  data  base  can  be  reduced,  this  translates  directly  into  savings  of 
time  and  money.  It  is  recommended  that  such  an  analysis  be  conducted 
for  the  applications  of  radar  image  simulation  most  often  used  by  the  U.S. 
Army  to  attain  its  mission  objectives. 

3.2.3  Data  Compression  Techniques 

Vast  amounts  of  data  must  be  p rocessed  for  a I 1 but  trivial  applications 
of  radar  image  simulation.  As  presently  structured,  data  bases  consist 
of  a point  in  a matrix,  at  least,  for  each  pixel  (picture  element)  in  the 
final  simulated  radar  image.  This  means  that  most  data  bases  for  opera- 
tional systems  are  exceptionally  l.irge  and  even  the  most  trivial  image 
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ItiincH  i lu)  is  i 11(^1  (I  i n.i  1 1' I V coiiipli'*,  Siinple  tt'ifuis  sdch  as  rotatio/is  of  data 

bast's  to  altfi  t lu‘  look  d i ri'C  t i on  (Mitihi  linci  are  t rcmendous  I y time 
cortsumiiui  and  expensive.  It  is  reromniended  that  both  lethni(|ues  tor  data 
compressit)n  and  alternate  merhods  toi  i n f oritia  t i on  storatie  and  retrieval 
l)e  i lives  t i 1)0 1 e<1 . Since  data  bases  n.i  I in  railai  im.uie  siimilation  have 
several  unique  tealiires  and  since  there  ore  several  critical  limitations  on 
them,  lechni()ues  o i’  d.iia  handlip.ii  and  compression  for  more  cieneral  imaqe 
proct'ssinu  are  not  necessaiily  optimal  and  may  not  be  applicable.  On 
the  other  haiul , proceiluies  lejected  tor  more  qeneral  imaqe  data  bases  may 
be  well  suited  when  aimed  at  r.idar  simulation  dal.i  bases  for  a specific 
application.  Thus,  i nves t i qa t i on  of  data  compression  techniques  which 
mi<|ht  be  viable  for  data  bases  for  radar  imaqe  simulations  should  be 
coupled  with  sensitivity  analyses. 

Alternate  mettiods  for  information  storaqe  and  retrieval  probably 
will  rt'quire  the  assembly  ot  a ‘.pecial  purfiosi-  com()uter  and  memory 
device  desiqned  strictly  for  imaqe  p r oces  s i nq  ap()  I i ca  t i ons  . The  (treat 
potential  value  of  radar  imaqe  simulation  as  a useful  tool  seems  to  argue 
that  this  i nves  t i (ia  t i on  needs  to  be  conducted. 

3.2.1)  Imaqe  Quality  Measurements 

The  c)uality  of  a simulated  radar  imoqe  or,  lor  th.it  matter,  a real 
iridar  imaqe  to  suit  its  intendr'd  purpose  must  presently  be  assessed  by  a 
human  be' i ik]  . This  de  1 1' rm  i na  t i on  ol  ipiality  is  very  subjective  and  relic's 
on  the  judqemi'nts  ol  diflerc'iit  peo|)lc’  who  have  diffc*rent  bases  of  expc’rience 
from  which  to  judge.  As  the  number  of  ap(i  I i ca t i ons  increases  for  both 
simulated  and  re<)  I radar  imaiies.so  also  increases  the  necessity  to  use 
more  objective  ir't.iqe  (lualitv  measurements.  Presumably,  these  objective 
itnarte  quality  mt'asuremc'n  t s would  u‘~e  a digital  computer  to  relate  the 
data  of  the  visual  reccird  (recordc'd  in  the  imarte)  to  application.  Cer- 
tain statistical  pr(>)ierfies  ot  trie  i rtraqt'  slmuld  be  used  to  predict  the 
usefulness  of  an  imaqe  to  satisfy  sin'cific  ob  i t'C  t i v('s  . Such  a set  of 
measurement  criteria  tio  trot  iircsi'ntlv  c'xist.  It  is  recommended  that  an 
investigation  he  conduett'd  to  dt'fine  such  a set  of  i nraqe  quality  measurement 
parameters . 

This  set  of  i macie  qtjolily  criteria  would  be  invaluable  to  the  recommended 
sensitivity  study.  In  fact,  a limited  study  of  the  type  recommended  here 
will  be  required  for  success  in  the  ''ensilivily  study  Moreover,  as  use 
of  radar  imaqe  simulation  increases,  necessity  will  decree  that  these 


23 


criteria  be  identiTied.  application  by  application.  The  opportunity  exists 
I o > OHM ) I i (1,1 1 r I lu'.c  i ml  i V i dii.i  I c I I ■ n I ■.  .iml  I n dc  I i iic  I hi'  i i i I cl  i .1  in 
advance  of  specilic  need.  Success  here  will  make  it  easy  to  relate  the 
value  of  radar  inuKie  simulation  to  each  new  application  and  will  allow 
dillerent  aiiencii's,  oi  labs,  or  pi'oph',  to  understand  how  a particular 
imacje  mav  satisfy  their  needs. 

3-2.5  Empirical  Data  Measurement  Program 

The  concept  of  radar  imaiie  simulation  has  been  shown  to  be  a valuable 

tool  for  a mu  1 1 i tilde  of  poten  t ia  I applications.  A major  obstacle  must  be 

overcome  before  the  promise  offered  by  radar  imaqe  simulation  can  be 

realized.  In  fact,  that  same  obstacle  is  a main  reason  why  radar  imagery 

is  not  used  to  the  fullest  extent  possible.  So  little  is  presently 

known  about  the  backscatter  properties  of  most  objects  in  the  terrestrial 

envelope  that  the  visual  record  of  this  parameter,  present  in  images  as  t 

j the  qreytone  variation,  is  not  properly  utilized  for  the  information  it  i 

t contains.  With  more  certain  knowledge  of  the  backscatter  properties  of  I 

1 . ' 

^ objects  in  a scene  would  come  a significant  increase  in  intelligence 

' gathered  from  each  image.  That  obstacle  is  the  necessity  to  have  backscatter 

data  available  for  as  much  of  the  terrestrial  envelope  at  as  many  microwave 

frequencies  and  polarizations  as  required  to  attain  the  mission  objectives 

of  the  U.S.  Army.  If  radar  imaqe  simulation  is  going  to  be  used  increas-  i 

inqly  by  the  U.S.  Army  (or  Department  of  Defense),  reliance  on  existing 

small  programs  and  the  mission  objectives  of  other  agencies  to  collect 

these  data  is  ill-founded.  It  will  require  many  years  and  good  luck 

for  this  need  for  backscatter  data  to  be  fulfilled  in  this  way. 

We  recommend  that  a five-pronged  study  be  implemented  to  gather  these 
required  data.  First,  the  major  projected  applications  for  radar 
imaqe  simulation  should  be  identified.  After  this  is  done,  the  specific 
microwave  reflectivity  categories,  frequencies,  and  polarizations  can 

' j 

be  identified  for  which  backscatter  data  will  be  required.  Second,  a | 

I I 

data  collection  system  should  be  developed  which  will  gather  these 
data.  Whether  this  system  is  a synthesis  of  existing  hardware,  or 
; designed  and  built  specifically  for  this  requirement,  or  both,  is  immater- 

ial. Very  probably  there  already  exists  hardware  in  the  defense  inventory 
which  could  be  used.  If  so,  this  good  fortune  would  reduce  the  cost  of 
this  phase  of  the  program  by  a substantial  amount.  Third,  a program 

2h 
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'.Iioiild  l)c  ill",  i .Hill  i IMP  1 1’liii'n  t I'll  1(1  (iillpct  the  lociiiirod  dilta.  Hio 

il.l  I .1  Ill  Ml  I ll  III'  I I 1 I I I'l  I I'll  till  .1'.  III.IIIS  . I '.I'.  I Ml.  I I V.l  I i .1  I i oil'.  Ill  I'.ll  ll  ( .1  I l'l|l  )l  V 
,1'.  .Ii'i.  I I'l'il  Illy  I 111'  i 11 1 .'lull'll  .ipp  I i I .1 1 i nil'.  . Ni"il,  ,1  piiiiiipiii  imisl  ho  ilosiiiiicii 
til  pi  iH  i"i'.  till'  il.it.i,  rills  p ii'ii  I ,1111  will  h.ivp  ti)  Drq.Hii/o,  CO  t a 1 oi)ue  , .iiui 
lilt'  tlH'Sf  il.it.i  ill  iiilil  i I i OP  ti'  proci'S'.  i lui  the  row  ilot.i.  Lost,  theorcticol 
I ec  1 1 oiiioiini' 1 i t .c.i  t t e r i ng  studies  iiei'd  Ici  he  iiiitioted  for  certoin  specific 
I eipi  i I eiiien  t s . It  i ■.  not  reosoii.tli  I e to  expect  to  meosure  bockscolter  doto 
lor  oil  pe  rmu  t o t i oils  of  the  'oriohles,  There  ore  too  man'/  comb  i no  t i ons  . 

This  is  where  theoretical  studies  show  their  value.  Theories  con  be 
(levelo(ied  to  extend  ond  extropolale  the  measured  data  for  the  cases 
which  ore  not  nieosurcd. 

This  proiiroMi,  or  some  other  realization  of  it,  is  required  if  the 
U.S.  Army  (or  Defense  Deportment)  is  to  have  on  active  <ind  viable  radar 
image  simulation  program.  It  is  not  the  specific  steps  of  the  program 
outlined  that  is  impoitonl,  it  is  the  collection  of  tfiese  data.  Without 
these  doto.  neither  can  the  concept  of  radar  image  simulation  nor  can 
radar  dot.)  itself  he  exploited  to  the  fullest  extent  possible  to  help 
t fie  Army  attain  its  goals  and  mission  olijec  fives. 

3.?.  6 Thcore  t i co  I .Models 

The  increasing  apti  1 i ca t i ons  for  radar  image  simulation  as  well  as 
interpretation  of  radar  imagery  require  an  even  larger  store  of  back- 
'.(Otter  data.  All  radar  i m.ige  simulations  must  use  some  model  for  the 
reflectivity  (iropeilies  ( hacksca t t cr ) of  the  objects  in  the  scene.  These 
d.tta  are  required  tiy  the  simulation  model  to  produce  ttie  greyscale  data  in 
the  simulated  image.  It  is  not  leason.thle  to  expect  to  measure  and  record 
the  hackscatter  data  lot  all  (lossihle  iiermu  t a t i ons  and  combinations  of  the 
variables;  frequency,  po I a r i /a t i oti , categories,  seasonal  changes,  and 
etc.  Theoretical  models  must  he  developed  to  extend  and  extrapolate  I he 
measured  ilata  to  cases  which  have  not  been  or  cannot  he  measured.  This 
is  a real  need,  not  a wli  i ms  i i .i  I musing.  Tfiese  theoretical  models  will  form 
an  ititegral  part  of  radar  imagt'  simulation  as  applied  tci  the  various 
a(ip  1 i ca  t i ons . 

3.?.  7 Other  Senscir  Sys  t ems 

As  it  is  helieved  that  the  trend  in  intelligence  gathering  systems  is 
I o use  e lec  f rn-o|i  f I ca  I systems  wfiicli  .ire  an  aggregate  of  a number  of 
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different  sensors.  It  is  retonmiended  Ihot  the  concept  of  simulation  be 
extended  to  these  other  sensors.  Radar  is  just  one  sensor  of  this  aqqre- 
qate.  Another  sensor,  for  example,  which  seems  to  have  great  potential 
value  is  the  FL I R (Forward- Look i nq  Infra-Red).  Other  sensors  also  exist. 
These  sensors  need  to  be  imideled  and  the  under  lyinq  phenomena  investigated. 
It  would  appear  that  transfei  func  lions  analogous  to  the  famous  radar 
equation  relating  the  various  aspects  of  sensor,  scene,  and  medium  can 
be  obtained.  Since  applications  toi  these  sensors  are  expected  to  increase, 
it  is  recof)»iiended  that  these  i nves  I i ga  t i ons  be  conducted. 

3.2.8  Develop  Applications  for  Radar  Image  Simulation 

The  general  concept  of  radar  image  simulation  has  been  developed.  What 
remains  is  to  extend  the  general  concept  to  specific  applications.  With 
development  of  each  application  will  come  unique,  new,  problems  which 
must  be  solved.  These  problems  will  run  the  gamut  from  feature  extraction 
for  building  data  bases  to  image  handling  problems.  Some  of  the  more 
important  problems  that  need  to  be  tackled  are  those  related  to  making 
operational  or  nearly  operational  systems  and  research  more  cost  effec- 
tive. The  studies  required  to  attain  this  goal  will  carry  over  into  the 
research  area  and  will  open  the  drjor  to  increasing  the  number  of  applica- 
tions for  radar  image  simulation. 

It  is  recommended  that  the  getieral  concept  be  developed  for  several 
specific  applications,  the  problems  identified  and  solved,  the  value 
determined.  Doing  this  will  document  the  value  of  radar  image  simulation 
and  will  increase  the  demand  for  solution  of  the  problems  presented  in  this 
section.  In  this  way  the  use  of  radar  image  simulation  can  be  accelerated 
and  exploited  to  the  fullest  extent  possible  to  help  the  Army  attain  its 
ob j ec  t i ves . 

3.2.9  Cultural  Scene  Simulation 

The  simulat ion  of  the  radar  return  from  cultural  objects  by  digital 

techniques  is  very  difficult.  Extensive  knowledge  of  the  geometric  and 

dielectric  properties  of  the  objects  is  required.  Also,  it  takes  a lot 

of  computer  time  to  use  these  data.  Alternate  techniques  for  the  sim- 

/ 

ulation  of  the  radar  return  from  cultural  objects  need  to  be  investigated. 


Ir 
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A prontisirni  pnlential  candiilate  is  offered  by  an  optical  filterinq 
technique'^.  This  optical  filterinci  model  needs  to  be  investigated  to 
determine  its  rantte  ol  applicability  together  with  alternate  processing 
teihnigues  to  rea  I i /i‘  the  modc'l.  For  instance,  the  model  miciht  be  reali/ted 
on  t h('  digital  computer  insti'ad  ol  the  optics  bench  by  dic|ital  bandpass 
filtering  techniques.  Techniques  need  to  be  devised  to  combine  the  optical 
radar  simulation  with  the  digital  simulation.  The  resulting  hybrid  radar 
image  simulation  model  would  combine  the  optical  simulation  for  cultural 
objects  with  the  diciital  simulation  for  distributed  objects.  Great  poten- 
tial savings  in  time  and  resources  are  offered  by  this  hybrid. 

In  addition,  other  techniriues  nec'd  to  be  investigated.  Development 
of  a number  of  techniciues  for  simulation  of  a particular  class  of  object 
would  sitnplify  the  simulation  process  and  reduce  costs. 

3.2.10  Validation  of  the  General  Radar  Image  Simulation  Model 


It  is  recommended  that  a study  be  conducted  to  produce  simulated 
imagery  of  a specific  area  for  comparison  with  real  imagery  of  the  same 
area.  The  area  selected  should  consist  of  as  many  of  the  variables 
encountered  in  radar  imagery  as  is  possible.  A good  data  base  should  be 
produced  for  this  area.  This  data  base  should  be  used  as  input  to  the 
simulation  model  and  the  output  simulated  radar  imagery  should  be  compared 
to  the  real  imagery.  The  results  of  this  study  would  establish  the  true 
potential  value  of  radar  image  simulation.  This,  then,  is  a way  in 
which  the  model  can  be  validated  and  will  establish  a baseline  of  quality 
for  radar  image  simulation. 


Frost,  V.  S.,  J.  L.  Abbott.  V.  H.  Kaupp,  and  J.  C.  Holtzman,  "A  Mathematical 
Model  for  Ter ra i n - Imag  i ng  Rada r and  Its  Potential  Application  to  Radar 
Image  Simulation,"  TR319-6,  Remote  Sensing  Laboratory,  The  University  of 
Kansas,  November,  1976- 
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ABBREVIATIONS  AND  ACRONYMS 


8 Angstroiirs  (10  meters) 

CRT  Cathode  Ray  Tube 
FLIR  Forward-Look i nq  Infra-Red 
PI  Principal  Investigator 
PP I Plan-Position  Indicator 
Pixel  Picture  Element 

RSL  Remote  Sensing  Laboratory,  University  of  Kansas,  Lawrence,  Kansas 

o°  Scattering  Coefficient  (Backscatter) . The  percentage  of  energy 
reradiated  in  a specific  direction  per  unit  area  by  a particular 
object. 

SLAR  Side-Looking  Airborne  Radar 

f>  Angle  of  Incidence.  Angle  between  antenna  electrical  boresight  and 
a local  vertical  from  the  antenna. 
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APPENDIX  A 


DIGITAL  MODEL  FOR  RADAR  IMAGE 
SIMULATION  AND  RESULTS 


The  following  technical  report  (TR  3I9“8) 
prepared  by  the  Remote  Sensing  Laboratory, 
The  Center  for  Research,  Inc.,  Uni'ersity 
of  Kansas,  is  included  in  this  report  to 
provide  the  rationale  for  the  discussion 
of  Sec  t i on  2 . 
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CRES 


-REMOTE  SENSING  LABORATORY 


ABSTRACT 


*: 

The  theory  and  implementation  of  a digital  c I osed -sys tern , radar  image  ■ 

simulation  model  are  reported.  Visual  results  of  several  SLAR  (Side- 
Looking  Airborne  Radar)  simulations  with  a slant  range  mode  are  presented. 

The  model  and  computer  software,  which  were  developed  at  RenKtte  Sensing 

Laboratory,  have  the  capability  of  producina  SLAR  or  PPI  (Plan-Position  ] 

Indicator)  imagery.  A test  of  the  integrity  of  the  model  and  software 
implementation  was  conducted  with  a data  base  consisting  of  computer- 

J 

generated  geometric  objects.  This  allowed  the  calculable  shadow  and  layover  i 

I 

response  to  be  employed  for  validation  of  the  operation  of  the  SLAR  irradel.  ] 

Though  the  shapes  of  objects  for  the  test  site  were  artificially  generated,  i 

the  dielectric  properties  encompassed  in  the  scene  were  derived  from  empiri-  i 

cal  backscatter  data.  This  input  reflectivity  information  originated  from 

the  agricultural/soil  moisture  data  bank  available  at  the  Remote  Sensing  1 

Laboratory. 

Pictorial  examples  reveal  that  the  simulated  imagery  realistically 
models  radar  effects  seen  in  real  imagery;  for  example,  layover,  shadow  and  i 

fading.  Examination  of  layover  and  shadow  on  the  imagery  produced  for 
various  flight  tracks  reveals  that  Che  simulation  package  is  performing 
well.  As  an  example  of  the  app 1 i ca t i on  of  radar  image  simulation  to  "real 
world"  problems,  a simple  study  of  the  visual  effects  caused  by  changing 
frequency  or  polarization  or  both  is  included.  The  frequency-polarization 
tests  indicate  the  applicability  of  simulation  for  optimum  discrimination 
studies,  feature  enhancement  tasks,  and  general  radar  system  design. 


1.0  INTRODUCTION 


A thporetical  model  has  been  developed  at  the  Remote  Sensinq  Labora- 
tory  in  an  effort  to  simulate  radar  imagery.’  The  result  is  a closed- 
system  model;  that  is,  a mathematical  formulation  which  encompasses  the 
phenomena  that  affect  the  known,  t ransmi t ted ,pu I sed  energy  and  subse- 
quently the  final  image  product  of  a radar  system.  The  implementation  of 
the  model  with  suitable  input  data  and  software  routines  for  a digital 
computer  has  been  accomplished.  SLAR  and  PPI  formatted  imagery  have  since 
been  successfully  simulated.  A description  of  the  radar  image  simulation 
theory  (in  particular  for  the  SLAR  model),  applications  of  the  model,  and 
experimental  results  are  included  in  this  document  in  addition  to  a copy 
of  the  digital  simulation  program  and  a brief  outline  of  the  necessary 
input  data. 

Radar  intage  simulation  (RIS)  techniques  have  been  applied  to  military 
tasks,  for  example,  interpretation  training,  navigational  aids  and  guidance 
systems  for  unmanned  airborne  vehicles'.  The  suitability  of  RIS  for 
technologically  advanced  electronic  guidance  has  been  recognized  in 
light  of  recent  testing  of  tactical  and  strategic  SLCM's  (Sea-Launched 
Cruise  Missiles)  and  ALCM's  (Air-Launched  Cruise  Missiles)  in  the  United 
States.  The  resources  exist  at  RSL  to  develop  simulation  models  for  po- 
tential hybrid  guidance  schemes  employing,  for  instance,  mu  1 1 i -spec t ra I 
sensors.  The  capability  to  predict  sensor  (and  successful  mission)  per- 
formance must  have  an  impact  on  the  design  of  such  future  systems  as  well 
as,  of  course,  use  in  providing  the  reference  information. 

RIS  also  has  importance  as  an  instrurnent  for  research  in  the  areas 
of  electromagnetics  and  scattering  theory  and  it  provides  a method  of 
predicting  and  optimizing  system  performance  when  information  concerning 
the  microwave  response  of  terrestrial  scenes  is  available.  Although  volume 
scattering  theories  and  hacksc:atter  data  fur  vegetation  have  existed  for 

' For  a more  comprehensive  discussion  of  applications  of  radar  simulation 
see,  "Military-Oriented  Applications  of  Radar  Image  Simulation," 
(authors  - J.  C.  Holtzman,  V.  H.  Kaupp,  and  J.  L.  Abbott),  Remote 
Sensing  Laboratory,  University  of  Kansas,  TR  319-10.  See  also 
Scientific  American,  February,  1977,  "Cruise  Missiles."  (author  - 
Ts  i p i s ) . 


several  ye.ti'.,  it  i •.  (loijhiliil  tii.it  this  know  I ptlcic  has  been  previously 
assembled  (that  is,  to  visually  p i t'd  i i t a i ad.ii  response!  in  the  mantn't 
presented  in  this  report.  lot  exani()ie  it  would  be  very  ditficull  lor  one 
to  postulate  or  ntentally  envision  the  differences  in  radar  imaqes  (for  a 
particular  complex  ground  scene  in  which  the  frequency  of  the  radar  was 
changed)  without  the  aid  of  simulation  of  the  imagery.  Thus,  the  point  of 
the  f requency /po I ar i za t i on  study  is  that  the  microwave  response  of  a 
terrestrial  scene  varies  rapidly  with  frequency  and  that  optimum  performance 
of  a discrimination  radar  can  be  sought.  Ar.jther  valuable  application  of 
RIS  is  estimation  of  seasonal  variations  on  the  radar  imagery,  particularly 
for  radar  guidance  systems  in  which  the  resolution  and  dynamic  range  of 
qreytones  (image  density)  are  fixed.  For  example,  it  would  be  important 
to  know  before  flight  whether  the  presence  of  a heavy  snow  fall  over 
a target  area  would  mask  the  otherwise  dominant  characteristics  of  the 
scene . 

This  work  is  organized  into  three  subtopics,  (1)  SLAR  model  theory, 
(2)  the  mechanics  of  producing  radar  image  simulations,  and  (3)  results 
which  include  simulated  imagery  From  several  experimental  flight  situations. 
The  ground  scene  for  this  work  is  an  artificial  data  base,  in  that  the 
elevation  data  was  assigned  by  a computer  to  describe  three  dimens'onal 
objects  on  a plane.  The  purpose  of  employing  a deterministic  scene  was  to 
allow  the  resea rchers  to  ca 1 cu 1 ate  known  shadow  and  layover  response  to 
a SLAR  at  a known  altitude  and  distance  from  the  near  range  of  the  data 
base  for  verification  of  accuracy  of  the  SLAR  simulation  program.  Since 
radar  return  also  depends  on  the  backscatter  characteristics  of  the  ground, 
empirical  signta  zero  data  were  associated  with  specific  areas  within  the 
scene.  The  significance  of  this  may  not  be  clear  until  it  is  realized  that 
the  simulated  imagery  will  accurately  represent  the  dynamic  range  of  image 
density  on  the  film  that  a real  radar  would  produce  from  the 
same  scene.  What  has  been  attempted  here  has  been  to  provide  a fine 
resolution  data  base  which  doubles  as  a test  mechanism  for  SLAR  simulation. 
Some  of  the  results  may  be  surprising  to  those  persons  not  familiar  with 
mu  1 1 i -f requency/pol ar i za t i on  imagery;  never t he  1 ess.  the  imagery  included 
reinforces  the  idea  that  radar  image  simulation  or  similar  modeling 
niques  may  have  many  interesting  app I i ca t ions . 


tech- 


Before  proceedinq  to  the  description  of  the  mechanics  of  simulating 
SLAR  images,  tlie  theory  which  is  the  basis  of  the  closed  system  riwd^l 

V 

will  be  presented  to  illustrate  the  determination  of  qreytones  (imaged 
dens i ty ) . 


2 . 0 IMAGE  SIMUIAT  I ON  _rHI  OKY 

In  an  iinaciinq  SLAR,  a short  pulse  of  mirrowavc  enrrqy  is  transmitted 
into  space  from  an  atitenna  whose  hori‘si(|hl  is  orthoi|onal  to  the  fliqht 
track.  This  energy,  confined  within  the  antenna  beam  solid  angle  (and 
side  lobes)  strikes  the  ground.  A portion  is  re  radiated  in  the  direction  of 
the  antenna  (amount  governed  by  backscatter  characteristics  and  geometry 
of  the  scene)  and  is  detected  by  the  receiver.  The  video  signal  can  be 
processed  in  various  ways.  for  example,  it  can  be  recorded  directly  on 
signal  film,  it  can  be  converted  to  digital  data  and  stored  on  magnetic 
tape,  it  can  be  displayed  by  intensity  modulating  the  beam  of  a CRT 
(Cathode  Ray  Tube)  and  photographed,  or  processed  by  other,  more  elaborate  methods. 
The  final  image  (the  visual  iecord  of  the  radar  return  signal)  depends  on 
variations  in  the  relative  strengths  of  the  signal  returned  from  different 
parts  of  the  area  imaged  to  produce  contrasts,  edges,  and  the  range  of 
image  brightness  (grcytoncs).  The  two  primary  factors  determining  the 
strength  of  the  signal  observed  on  the  final  image  produced,  and  consequently 
the  brightness  of  an  image  point,  are  geometry  and  dielectric  properties. 

Simulation  of  radar  images  may  be  accomplished  by  mathematical 

operations  on  the  knov;n  parameters  of  the  radar  to  he  modeled  and  site 

to  be  imaged  to  form  a visual  display.  The  general  technique  utilized 

2 

m this  document  is  reported  by  Holtzman,  et.  al.  In  this  method  the 
radar  equation  is  used  to  relate  target  empirical  backscatter  coeffi- 
cients to  relative  image  greytones.  A general  form  of  the  greytone 
expression  is  developed  in  this  section.  Also  included  in  the  simulation 

model  are  (1)  the  effects  of  relief  and  tilting  of  resolution  cells,  (2) 

2 

layover  and  shadow,  and  (3)  radar  fading. 

The  general  model  for  mean  power  received  is  based  upon  the 
radar  equation  which  may  be  expressed  as^ 

2 

Holtzman,  J.  C..  V.  H.  Kaupp,  R.  L.  Martin,  E.  E.  Komp,  and  V.  S.  Frost, 

"Radar  Image  Simulation  Project:  Development  of  a General  Simulation 
Model  and  an  Interactive  Sin)ulation  Model,  and  Sample  Results," 

TR  23'*-13.  Remote  Sensing  Laboratory,  The  University  of  Kansas,  Feb.,  1976. 

^ Moore,  R.  K. , Remote  Sensing  Manual.  (Editor  - Reeves),  Chapter  9, 

American  Society  of  Photogrammci ry , 1975. 


(I) 


p 

r 


P^  ■ ■ ' A ■ r/  ( i) 

{'4'  ) 


where  P^  = Average  return  power  received  at  the  antenna  terminals; 

P^  = Average  transmitted  power; 

= Scattering  coefficient  of  the  ground  spot  resolution  cell 
(function  of  local  angle  of  incidence,  0 , the  incident 
wavelength  (‘1,  t ransm i t/rece i ve  po 1 ar i 7a t i on , and  surface 
parameters) ; 

0 = Radar  incidence  angle  to  surface; 

Gt&l- Antenna  gain  in  the  direction  of  AA  (function  of  the  radar 

system  mcjdeled  and  assumed  to  be  identical  in  the  return  direc- 
tion) ; 

,<  = Wavelength  of  the  incident  wave; 

R = Distance  to  the  ground  spot  resolution  cell  (A  function  of 
altitude  and  look  angle); 

AA  = Area  of  the  ground  spot  resolution  cell. 

The  area  'A  of  a reso 1 ut i on  ce I 1 for  a pulsed  radar  can  be  modeled 

2 

as 


where:  w = Size  of  resolution  cell  in  the  along-track  direction; 

= Size  of  resolution  cell  in  the  cross-track  direction; 
h'  = Height  difference  between  the  cell  and  the  radar; 
if  = Antenna  beam  width; 
g = Radar  incidence  angle; 

('  = Sional  pulse  width; 


0^  = Local  slope  of  resolution  cell  in  the  along-track; 
0^  = Local  slope  of  resolution  cell  in  the  ac ross- track . 


Holtzman,  J.  C.,  V.  H.  Kaupp,  R.  L.  Martin,  E.  E.  Komp , and  V.  S.  Frost, 
"Radar  Image  Simulation  Project.  Development  of  a General  Simulation 
Model  and  an  Interactive  Simulation  Model,  and  Sample  Results," 

TR  23^-13,  Remote  Sensing  Laboratory,  The  University  of  Kansas, 
February,  1976. 


If  the  averaqe  return  pov/er  P from  a particular  category  of 

^ I 

siatterer  ( f, ) at  a iiaitiiul.ii  aiuile  of  ipcidetice  ("i  ) for  a calibration  sys- 
tem is  known,  then  equation  (I)  can  be  rewritten; 

Pt  ...-P,,  I ■ '.A,  • r, 


Equation  (1)  can  be  used  to  define  the  averaqe  return  power  Pj. 

^2 

from  the  same  category  of  scatterei  (C)  at  a different  angle  of  incidence 
("2^  for  the  radar  system  being  modeled. 


Vc‘'’o’  ' '“a  ■ s’"’?'*; 


r - , V-’t 

C2 

Dividing  Equation  fit)  by  Equation  (3)  gives: 

P^.  n"  (fl.  )AA-G-^(0  )>, 

T.,  c sf-  2 2 2 2 1 

p _ p ^ £ (5) 

% ^'1 

Thus  the  averaqe  return  power  for  each  point  in  an  image  (for  any  parti- 
cular category  of  scatterer)  can  be  found  if  the  average  return  power  for 

one  point  belonging  to  a category  is  known  (0°^  absolute)  in  addition  to 
backscatteri ng  coefficients  and  angles  of  incidence  for  both  points. 

The  visual  presentation  medium  modeled  was  photographic  film.  An 
image  density  which  is  related  tc'  the  average  return  power  is  displayed 
on  real  SLAR  images.  This  image  density  is  often  called  a qreytone*, 
which  is  a relative  measure,  that  is,  it  was  produced  with  respect  to 
some  calibration  reference.  I \ a digitized  image,  the  greytone  repre- 
sents a specific  image  density  level  within  the  possible  dynamic  range 
on  the  image.  The  return  power  (related  to  video  intensity),  properties 
of  the  film  and  photographic  processing  methods  determine  the  photographic 

if 

density  (D)  on  the  film  by 


Goodman,  J.  W. , I n t roduc t i on  t o Four i er  Optics,  Chapter  7,  McGraw-Hill, 

1968. 


0 = f 


where:  y = Gaimiia  of  the  film; 


I,.  = A positive  constant  which  depends  upon  the  exposure  time  and 
on  the  film  processing  and  development; 

I = I maqe  i n tens i ty ; 

D = Photographic  density. 

This  relationship  holds  true  In  the  linear  portion  of  the  film  dynamic 
range.  If  a linear  radar  receiver  is  assumed,  the  the  intensity  (t)  is 
directly  related  to  the  average  return  power  (P^)- 

I = M • P (7) 

r 

where;  M = A proportionality  constant. 

Rewriting  (f>)  to  incorporate  (7)  gives 


D = >log,oP^  + log^pk  + flog,QM 


This  result  will  be  used  for  the  final  greytone  expression. 

4 

The  photographic  density  (0)  is  defined  by 


D = )og,o(-) 


where  the  film  transmittance  i is  given  by 


with:  1^  = Transmitted  intensity; 

I = Incident  intensity, 
o 

Rewriting  Equation  (3)  produces 


D = iog,o(-i^) 
t 


Goodman,  J,  W. , Introduction  to  Fourier  Optics,  Chapter  7.  McGraw-Hill, 

1968. 


The  Morn.<ili/ed  v.i  I ue  ot  t fie  density  (1>)  is  the  .|i‘’ylone  level.  1 h.' 
iiiMMO  simulations  were  d i q i t i /od  into  8-bit  binary  words  for  display, 
allowiiui  256  distinct  qreytones.  If  the  optical  density  used  in  a simu- 
lation is  in  the  ranqe  from  0 to  X,  the  lineal  portionof  the  film  dynamic 
I anqe  is  required  to  lie  of  sufficient  s i 7e  that 


In  the  simulated  imaqery  produced  at  RSL,  a qreytone  of  25b  was 
associated  with  white  and  zero  with  black.  This  is  the  expression  for  a 
positive.  This  siqnal  is  then  used  to  intensity  modulate  a CRT  (a  posi- 
tive imaqe).  This  i maqc  is  then  photoqraphed  (a  neqative).  The  neoative 
is  then  used  to  produce  a photo  (a  positive)  which  is  the  final  product. 
The  qreytone  level  G correspondi nq  to  a density  D was  qiven  by 

G = — ■ D (13) 


The  qcneral  relationship  between  the  average  return  power  and  qreytone 

level  for  each  resolution  cell  in  tlie  simulated  image  is  obtained  by 

rewritinq  Equation  (8),  i ncorpora t i nq  Equations  (5)  and  (13): 

r .r  (n  , )aA  G/(n  ) ,1  2r 

2 bb  \ T ^ ? ''7  22  1 I 

's ’ "q  * I’”"'” 


+ loq 


10  k 


This  equation  states  that  the^relative  qreytone  values  for  eacti  point 
in  an  imaite  corresponding  to  any  anqU'  ot  incidence  ol  a category  of  scatt- 
('ters  can  be  obtained  if  the  qrey'one  and  angle  of  incidence  is  known 
for  one  discrete  point  in  any  category.  This  equation  is  the  general  result 
which  establishes  the  relationship  between  the  average  return  power  and 
relative  qreytone  levels  for  a simulated  radar  image.  It  should  be  pointeil 
out  that  on  a radar  imaije  and.  thereby,  on  a simulated  radar  image,  the 
important  parameter  is  the  relative  qreytone  level  between  points  (in  the 
image  and  not  the  absolute  values  of  the  qreytone.  Since  this  is  true. 
Equation  (|4)  can  be  used  to  establish  the  relative  qreytone  levels  between 
all  categories  of  scatterers  included  in  an  image  provided  that  the  absolute 


where : 


V. tines  of  the  srolterint)  coe  f t i f i en  t , and  the  appropriate  anqies  of 

incidence  foi  all  points  are  known  and  provided  that,  in  addition,  the  nrey- 
tone  of  one  point  in  the  i maqe  is  known.  The  more  general  result  is,  then; 

r.  - h • ' Ion  ; I ] H'-) 

'*^2 

+ loq  + . loci  I 

10  kj  ' '"‘'lO^M,  ' 

. 

where:  = Greytone  level  of  some  point  helonqinq  to  category  'c' 

*"2  anfi  having  the  scattering  coefficient,  ‘'^2° 

at  the  appropriate  antjle  of  incidence,  by  a radar  sys- 

tem; 

Gj^  = Known  greytone  level  fora  particular  point  belonging  to 

^1  category  'a'  and  having  the  scattering  coefficient,  i'  , ° (B  , , ] 
measured  at  the  appropriate  angle  of  incidence  B jj  | . by  a 
cal i bra  ted  sys  t em. 

Equation  (15)  is  the  general  theoretical  result  which  has  been  implemented 
ina  computer  simulation  package.  Implicit  in  this  equation  are  the  effects 
of  geometry  which  cause  the  phenomt?na  of  shadow,  layover  and  reflection  of 
microwave  enerqy  (in  accordance  with  local  tilt  of  resolution  cells). 
Reference  [2]  contains  a ntore  complete  treatment  of  these  effects. 


2 . I Fad i n g 

Equation  (I5)  qives  the  impression  that  a deterministic  process  is 
occurring.  The  interaction  between  the  radar  and  a surface  is  a random 
process  to  some  extent,  and  the  mean  value  of  the  greytone  level  for  a 
given  set  of  parameters  is  described  by  Equation  (15). 

If  one  assumes  that  the  process  can  be  modelled  as  additive,  i.e. 
signal  = mean  + noise,  a model  for  the  true  statistical  properties  of  the 
phenomena  can  be  derived.  The  return  power  froni  a single  scatterer  has 
a noise-like  characteristic  that  follows  a Rayleigh  distribution^.  This 


Moore,  R.  K. , Radar  for  Geoscience  Instrumentation.  Chapter  5.6,  Geoscience 
Instrumentat  ion,  Ce  . A . Wolfes,  Ed  i tor)",  John  Wl  ley  and  Sons. 


distribution  is  characterized  by  equal  mean  and  standard  deviation.  The 
probability  for  N‘  such  scatterers  follows  a chi-squared  distribution  with 
2N  ' deqrees  of  freedom.  The  mean  and  variance  of  chi-squared  distribu- 
tion are  given  by : 

mean  = deqrees  of  freedom  = 2N'=  p 
variance  = 2 (mean)  = ' 


standard  deviation  = 


mean 


For  radar,  the  number  of  independent  saniples  in  a given  resolution 
cell  is  given  by 


^ = ^h'j.  ^ 
L )cosB 


(16) 


where:  0 = Radar  incidence  angle; 
w = Azimuth  resolution: 

L = Horizontal  aperture  length  of  antenna  (L  - — ) ; 

(J* 

h'  = Effective  altitude; 
i|>  = Antenna  beaniwidth; 

> = Wavelength  of  microwave  energy. 

Therefore,  the  standard  deviation  due  to  fading  can  be  represented  by: 
mean  _ -g-  _ _LL  (1 7 ) 

v'N  ^ ►¥ 

For  a particular  cell,  the  return  power  (P^)  can  be  represented  as: 

P = II  + I • RN  (18) 

r X 

v-jhere:  RN  is  a random  variable  which  has  the  d^i  r«d  ■ ha  racte  r i s t i cs  (chi- 

square  distribution  with  mean  and  variance  described  above). 

For  computational  purposes  the  random  variable  (RN)  can  be  normalized 
to  zero  mean  and  variance  of  one  leading  to; 

P = i:  + X RN  = i.  (1  + — ) (19) 

’’  v'N  v'N 

^ Moore,  R.  K. , Radar  For  Geoscience  Instrumentation,  Chapter  5.6,  Geoscience 
Instrumentat ion,  (E.  A.  Wolfes,  Editor),  John  Wiley  and  Sons. 
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iiii.mo  inlonsiiv  .i  spoi  i I i od  SI  AK  would  piddiico  cm  rospond  i mi  to  <1  piii  * 
ticul.ir  tniqol  . 110.1.  tiooiiiotiic.il  otiocis  I h.i  I I'liti'i  t ho  ilo  lorini  n.i  t ion  of 
t.iiiiot  .110,1  (rosoliition  coll  i /o  .ind  loc.il  o i i on  t .1 1 i on ) .no  tro.itod 
I iiioiously  in  tlio  diiiit.il  s i imi  I .i  t i on  iiiodol  .ind  .no  t tioromih  1 y covorod  in 
Rotoii'nco  I2I.  Till-  SLAR  s i iiui  1 .1 1 i on  pioiii.iiiis  .110  .ihio  to  t1c.1t  hoth  the 
c.iso  ol  littio  or  no  .ind  ol  •,  i qii  i I i c.in  t roliol.  Socond'ordi'i  olt'octs 
such  .IS  imiltip.ith  which  in.iy  occni  diio  to  loc.il  iioomo  I 1 y li,ivo  not  yot 
hoon  i nip  1 onion  t od  in  t ho  SIAR  proni.iiii. 

The  succoodiiui  sections  di'.cnss  the  iii.icrollow  ch.irt  ol  t ho  SLAR 
pi  oqi  .nils , I ho  input  d.it.i  which  ni.iy  ho  .idjiistod  hy  I ho  user  .iccoril  i ni|  to 
I ho  .ipp  I i C.I  t I on  , .iiid  the  d.i  t .1  h.iso  th.it  h.is  hoon  ,111  oxpor  i nion  t .i  I tool  tor 
V.1 1 i d.i  I i on  ot  the  c losi'd  systoii!  niodo  I tor  d i q i t ,i  1 s i mu  1 q t i on  ol  r.id.ii 
ini.iqos.  A copy  ol  the  SLAR  proqr.inis  is  included  <is  Appondix  I. 


Ilo  I I /ni.in , ,1.  I'..,  V . 11.  K.nipp,  R.  1.  M. lit  in,  L.  L.  Konip , ,ind  V.  S. 

Frost,  "R.id.ii  lin.iqo  S i mu  1 .1 1 i on  I'lojocl:  Dovolopiin'iit  ol  ,1  rionor.il 
S i mu  1 .1 1 i on  Model,  .ind  S.imple  Results,"  TR  PVt-l},  Romoti'  Sonsinq 
L.ihor.i  t ory  , The  tin  i vi' 1 i I y ol  K. ins. is,  Lohiii.iiy,  I'V/h. 
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3 . 0 SLAR  S LMU  L_AT  ion  P ROC  E S_S 

Tho  SLAR  simulation  packaqe  (a  copy  is  included  as  Appendix  I) 
consists  of  four  separate  programs  and  associated  subroutines  as  given 
in  Figure  1.  They  can  be  run  separately  so  that  (1)  the  failure  of  one 
module  does  not  affect  successful  comoletion  of  the  previous  programs,  and 
(2)  experimentation  can  be  conducted  in  one  program  for  a particular  radar 
being  modeled,  thus  avoiding  the  cost  of  rerunning  the  entire  package. 

For  example,  it  might  be  desired  to  change  a reference  greytone  or  the 
dynamic  range  in  the  simulated  image.  This  sturcture  for  the  software 
would  allow  SLAR  Geometry,  SLAR  Slope,  and  SLAR  Shadow  (the  first  three 
programs)  to  be  run  once,  and  the  output  tape  to  be  used  for  several 
successive  experimentations  of  SLAR  Greytone  (the  fourth  program). 

MACRO-FLOW  CHART  OF  SLAR  SIMULATION  PROCESS 
( START  ) 


Compute  Rcsclution  Cells 
(SLAR  Geometry) 


L 

~ — — ^ - 


i.l  Input Puic'iniut  e-is 

Foul  .irinys.  HADAK,  aROUNO,  KLFLR  .irul  PROCES  contjin  locutions  lor 
il.itu  which  vury  with  the  desired  dpi)  1 i cn  t i on  of  the  SLAR  simulation  and 
may  he  adjusted  hy  the  user.  Descriptions  of  the  arrays  are  contained  in 
Fiiiure  2 wtiich  is  a summ.iry  ol  all  (larameters  (summarized  from  propram 
listinu  in  A()pendix  I)  which  are  chanpeable.  The  asterisks  denote  excep- 
tior.s  as  explained  at  tiu'  conclusion  of  tfie  listinq.  The  user  must  have 
a minimal  knowledqe  of  the  interdependence  of  the  radar  and  qround  para- 
meters before  attemptinq  to  alter  the  values  in  these  matrices.  For 
example,  a clianqe  in  the  parameter  RADAR(6)  may  necessitate  compensation  in 
GROUND  (1),  (2).  (5).  (10)  and  RADAR  (7),  (8).  and  (9)  parameters. 

5 . 2 Da^.i^ 

An  artilicial  data  base  (scene  consist imi  of  qeometrical  solid 

shapes)  was  used  to  test  the  operation  of  both  the  simulation  model  and 

the  software.  A data  base  of  qec^metrical  sha(>es  was  used  for  these 

IMirposes  because  of  the  known  i nfiut /output  response  characteristics  of  sucfi 

a data  fiase.  Ttie  input  was  known  exactly  and  the  output  was  deterministic 
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and  calculable.  The  artificial  data  base  described  by  Komp  is  a composite 
ol  Ihrei'  dimension.il  ohjiulv  Tln'  ma  l i i x of  elev,)lii>n  ,ind  caleqoiy  d.il.i 
( 1 1 u I 'd  III  ,1  I I iiiipii  1 1 ' I wi  u 1 1 1 1 11  I n 1 1 pi  r.  I I 1 1 >1 1 ) wh  I I h I i > i ms  I I n ■ d . 1 1 . i h.  i .i  ■ 

is  700  by  1000  points  in  extent,  represent inq  20  by  20  foot  resolution. 

Thus,  the  data  base  appears  to  be  Ut.OOO  by  20,000  feet.  Known  empiri- 
cal backscatter  characteristics,  in  the  form  of  s i qma  zero  versus  anqie  of 
incidence,  for  ten  separate  veqetation  cateqories,  have  been  assiqned  to  the 
qround  scene.  The  value  of  such  a data  base  with  completely  specified  die- 
lectric and  qeometric  properties  is  that  it  provides  us  with  a means  to  vali- 
date the  SLAR  simulation.  The  calculable  response  (e<i.,  extent  of  sliadow  and 
l.iyovi'r  due  to  each  objc'ct)  of  the'  qeometric  solids  to  various  fliqtit 
tracks,  freiiuency/tiolar  i zat  ion  combinations,  near  and  far  ranqe  depression 
anules,  etc.,  allows  the  final  imaqe  product  to  he  judcied  objectively. 
Naturally,  a hi<|h  resolution  (or  larqi'  matrix)  dat.i  hast'  can  he  qenciated 
I'y  a computet;  liowever  , the'  numf'c'r  of  data  points  used  will  be  limited  by 
t lie  specific  aiiplication  and  computer  ri'sources. 

^Komp,  F.  D.  , V.  If.  Kaupt',  and  .I.  F.  Moltzman,  "Construction  of  a Geonu'tric 
Data  Rase  lor  R.idar  Imaqe  Simul.ilion  Studies,"  TR  3I9-I,  Ri'iixite 
Sensinq  laboratory,  Ffie  University  of  K.tnsas,  .)uly,  1976. 


The  d.Ttn  h.isi-  ot  qeiimetric  solirl^  is  shown  in  Ficiuie  3.  The  flat 
ctieckerbOiird  roncep  r no  I I y represents  oreas  of  several  vegetation  types. 

The  objects  consist  ot  lu'ini  stitiei  es  . i v I i nde  is  . pyramids  and  rectaiuiulai 
para  1 le lop i peds , with  heiqhts  and  dimensions  drawn  to  scale.  The  variety 
of  shapes  present  many  different  local  tilts  to  the  radar  and  will  cause 
qreytone  cbanqes  at  the  sidt's  of  objects.  By  makinq  the  data  base 
large  and  with  resolution  much  finer  than  the  siimulated  radar  image,  it 
was  possible  to  study  many  properties  of  the  simulation  software  (e.q.,lhe 
formation  of  resolution  cells  by  the  SLAR  package  and  the  layover  effects). 

Figure  h shows  the  relative  position  and  elevation  of  the  many 
objects  on  the  data  base.  Thi;  maximum  variation  of  relief  is  2000  feet, 
therefore,  a wide  range  in  the  aiiwunts  of  shadow  and  layover  are  expected 
to  be  evident  in  simulated  imagery  of  this  scene. 

Figure  5 illustrates  the  microwave  reflectivity  category  assign- 
ments for  the  objects.  Comparing  Figures  4 and  5,  it  is  seen  that  the 
same  vegetation  covers  several  objects  of  differing  height.  The  value  of 
empirical  backseat  ter  data  for  simulation  is  readily  apparent,  for  if  it 
were  not  available,  complicated  analytical  expressions  derived  from  volume 
scattering  theories  would  have  been  necessary  as  input  data.  Backscatter 
data  at  increments  of  one  degree  between  zero  and  ninety  degrees  (incidence 
angle)  were  stored  (for  each  vegetation  type)  on  input  file  tapes. 

The  SLAR  simulation  produces  two  types  of  output:  (1)  a "shadow 
map"  (hard-copy  output  of  the  SLAR  package)  and  (2)  a magnetic  tape  compat- 
ible with  a digital  computer  controlling  a display  medium.  The  shadow 
map  consists  of  a large  matrix  of  characters,  zeroes  representing  shadow 
areas,  nines  indicating  no  data,  and  all  other  numbers  the  sum  of  the  data 
points  falling  within  a reso 1 u t ion  ce  I I bin.  Both  types  of  output  are 
illustrated  in  the  following  sections  for  a radar  system  operating  at  8.6 
GHz  with  HH  polarization  (with  the  exre(>tion  of  the  f retiuency/po  1 a r i za  t i on 
study  results).  Thus,  the  shadow  map  indicates  the  shadow  as  well  as 
layover,  and  sums  the  number  of  I ransposed  cells. 


4.0  DIGITAL  SLAR  SIMULATION  RESULTS 


The  dato  has<*  information  (rJata  point  location  in  the  array,  hciqht 
and  catecioiy)  crrntained  on  niaiinetic  tape  ^ qenerated  by  the  methcKl  of 
Reference  [6]  was  used  as  the  qround  scene  for  the  SLAR  simulation  programs 
Two  sets  of  visual  displays  were  produced,  shadow  maps  and  radar  image 
s imulat ions. 

The  map  shown  in  Figure  6a  illustrates  the  effects  of  layover  and 
shadow  that  would  occur  for  a SLAR  operating  at  8.6  GHz  HH  flown  at  a 
40,000  foot  altitude  and  20°  near  range  depression  angle.  Comparing  this 
map  to  Figures  3 and  4 it  is  seen  that  the  level  checkerboard  pattern 
lies  in  the  upper  portion  of  Figure  6a.  The  symbols  on  the  map  indicate 

the  sum  of  the  resolution  cells  whr)se  radar  returns  simultaneously  reach 
the  SLAR  receiver.  Zeroes  represent  shadow,  or  lack  of  appreciable  return, 
and  nines  are  employed  to  show  no  data.  Shadow  and  layover  effects  do 
not  occur  on  the  flat  terrain,  just  as  would  be  predicted.  Careful  study 
of  Figures  3,  4,  and  6a  reveal  that  the  SLAR  package  is  accurate  in 

handling  radar  geometrical  phenomena.  This  particular  version  of  the  sim- 
ulation limited  the  number  of  resolution  cells  laying  over  into  one  bin 
to  a sum  of  six.  Examination  of  the  map  shows  that  the  maximum  number  of 
points  to  enter  one  bin  was  five,  thus  the  above  restriction  did  not  limit 
the  performance  of  the  programs. 

Figure  6b  simulates  the  shadow  and  layover  effects  for  a SLAR 
(8.6  GHz  HH)  at  40,000  feet  and  a near  range  depression  angle  of  40°. 
Shadows  are  shorter  and  layover  more  pronounced  for  this  case  than  for  the 
previous  one.  The  centrally  located  pyramid  in  6b  shows  a larger  area  of 
layover.  However  each  bin  contains  a maximum  of  three  resolution  cell 
radar  returns.  The  smaller  area  of  layover  in  Figure  6a  has  bins  with 
up  to  five  resolution  cell  returns  summed.  The  more  distinct  outlines  of 
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Martin,  R.  L.  Martin's  Thesis,  "SLAR  Simulation  and  Applications,” 

University  of  Kansas,  l‘176.  (This  document  addresses  problems  of 
handling  in)>ut  and  out  put  data,  tape  compatibilities  of  machines, 
and  system  rou tint's.) 
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Komp , E.  0.,  V.  H.  Kaupp,  and  J.  C.  Holtzman,  "Construction  of  a Geometi ic 
Data  Base  for  Radar  lma<|e  Sinitjiation  Sti/dies,”  TR  319*2,  Remote 
Sensing  Laboratory,  The  University  of  Kansas,  July,  1976. 


6b  allow  each  object  to  be  iilentified  in  conjunction  with  Figures  3 
and  . 

It  can  be  seen  that  radai  interpretation  for  niiliiarv  personnel 
could  bt'  facilitated  by  the  use  of  shadow- 1 ayover  maps  similar  to 
those  presented.  The  shadow  map  allows  one  to  understand  differences 
in  image  formation  between  radar  and  conventional  photographic  processes. 

The  distinctions  arise  because  the  radar  emits  pulsed  energy  such  that 
the  return  energy  can  be  plotted  versus  time,  the  lack  of  return  power 
causing  shadow,  and  the  propagation  time  determining  the  placement  of 
relative  greytones  due  to  return  signals.  However,  in  the  case  of  optical 
photographic  imaging,  the  source  of  illumination  (the  sun)  emits  energy 
continuously,  the  reflected  energy  giving  rise  to  a video  signal  which  cannot 
he  charted  versus  time  without  considerable  ambiguity  as  to  the  time  ori- 
gin of  the  transmitted  signal.  Thus  the  continuous  superposition  of 
reflected  energy  masks  shadenv  formation  (in  the  sense  of  lack  of  return 
signal)  when  the  photographic  imaging  platform  is  in  the  path  between 
the  illuminator  and  the  object  to  be  sensed.  The  handicap  imposed  by 
familiarity  solely  with  photographic  images  is  difficult  to  overcome, 
and  this  can  only  be  accomplished  through  practice.  Experience  with 
shadow  and  layover  for  various  flight  situations  simulated  with  the  aid 
of  an  artificial  data  base  could  be  extremely  useful  for  overcoming 
the  urge  to  think  in  terms  of  photographic  phenomena.  The  shadow  map, 
by  the  very  fact  that  it  contains  a limited  amount  of  guantified  infor- 
mation about  the  processes  occurring  in  radar  terrain  imaging,  provides 
the  necessary  elements  for  (jnderstanding  shadow  and  layover.  Despite 
the  fact  that  the  geometrical  data  base  is  not  representative  of  common 
landscapes,  the  corresponding  maps  (jenerated  by  simulation  are  valuable. 

Simulated  radar  images  for  a SLAR  altitude  of  4000  feet  and  a near 
rantie  depression  angle  of  40  degrees  (far  range  ten  degree  depression  angle) 
were  produced  for  a f reguency/po I a r i 7a  I i on  study.  The  accentuation  of  sha- 
dow and  layover  was  caused  by  the  presence  of  objects  on  the  data  base 
up  to  2000  feet  high  and  the  extrentely  low  flight  altitude.  The  results 
are  shown  in  Figure  7 for  six  different  f reguency/po I a r i za t i on  combinations. 
Several  interesting  radar  effects  can  be  spotted  by  an  untrained  eye:  (I) 


near-ranqe  comp icss ion ; (2)  backseat  lor  cateqory  d i scr I mi na t i on  and  vari-  j 

ability;  (3)  samc-ca teqory  blocks  of  the  checkerboard  chanqinq  qreytene  ] 

across  the  map  because  of  clianqes  in  s i qma  zero  data  with  (>  roqress  i n(|  ] 

an<tle  of  incidence',  ttieta;  (h)  selective  hi  i uh  1 en  i nq  ol  objects  with  polaii-  I 

zation  adjustment;  (5>)  obviously,  layover  and  shadow;  and  (6)  loss  of  siqnal  1 

due  to  local  lilt  ot  resolution  cell  on  sides  of  cylinders. 

It  is  possihli*  that  a fine  resolution  radar  (and  correspond i nq I y 
liner  data  base  resolution)  would  indicate  the  presence  of  hard  taiujets 

on  the  qeometrical  data  base  due  to  some  of  the  corners  and  edqes  c>f  the  ' 


shapes  involved,  althouqh  the  catectories  strictly  represent  vonetation. 
However,  this  would  not  be  likely  if  the  real  SLAR  produced  60  foot  resolution 
as  simulated,  for  the  random  averaqinq  processes  across  a resolution  cell 
would  somewhat  mask  the  effects  of  return  power  from  the  hard  tarqets 
(eq.,  top  edqes  of  cylinders  and  parallelepipeds).  Treatment  of  this  data 
base  (with  computer  generated  resolution  as  fine  as  desired)  as  a composite 
r>f  distributed  tarqets  is  accurate  as  lonq  as  coarse  (eq.,  6<)  x 6o  foot) 
resolution  is  employed  in  the  simulated  imagery. 

Assume  that  the  flight  mission  of  a real  SLAR  v^;as  to  produce  cateqory 
information  over  flat  terrain  and  that  a low  altitude  constrained  the  near 
and  far  range  depression  angles  such  that  only  the  middle  third  of  the 
images  in  Figure  7 were  sensed.  Such  a mission  would  be  fruitful  if  the 
particular  radar  system  was  operating  at  8.6  GHz  with  HH  polarization. 

With  any  other  combination  of  f requency/pol ar i za t i on  shown,  the  information 
gathered  would  not  justify  the  expense.  Radar  parameters  have  already 
been  studied  in  this  manner  for  the  Earth  Resources  Shuttle  Imaging  Radar 
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to  predetermine  discrimination  ability  . Whether  the  imaging  system  is  a 
SAR  or  SLAR,  waste  and  inefficiency  in  related  missions  can  be  avoided  by 
forethought;  a helpful  aid  would  be  the  simulation  of  radar  images  from  data 
base  with  relief  and  category  (hopefully,  empirical  backscatter)  information, 
be  it  a real  site  or  the  construct  of  a geographer's  imagination. 

Rather  than  simiilatinq  flight  at  4000  feet  in  the  midst  of  terrain 
varying  in  height  from  0-2000  feet,  an  altitude  of  40,000  feet  might  be 
more  reasonable  for  a SLAR.  Experime-nts  have  been  conducted  with  two 


Bush,  T.  F.,  "Cropland  Inventories  Using  a Satellite  Altitude  Imaging 

Radar,"  Ph.D.  Dissertation  to  t)e  published  (Spring  1977),  University 
of  Kansas. 
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Mi<|ht  qeonietries  as  shown  in  Figure  8.  The  extent  of  shadow  and  lay- 
over will  differ  in  the  cases  of  20  and  2k  degree  near  range  depression 
angles.  Figure  9 illustrates  that  shadiiws  will  be  longer  in  the  20“ 
example,  and  that  data  points  layover  greater  distances  toward  the  flight 
track  in  the  2k°  depression  angle  case,  especially  in  the  near  range  (this 
assumes  ground  range  presentation).  However,  more  brightening  will  occur 
in  the  layover  areas  of  the  20"  case  due  to  the  greater  density  of  cells 
entering  bins  on  the  near  range  side  of  tall  objects. 

The  resulting  simulated  radar  images  are  shown  in  Figures  10  and 
II.  Close  examination  of  lengths  of  shadows  and  layover  are  aided  by 
employing  the  checkerboard  as  a distance  marker.  Allowing  for  the  dis- 
tortion in  photographs  of  the  IDECS  visual  display,  the  shadows  are 
longer  in  Figure  lo  and  data  points  lay  over  greater  distances  toward  the 
flight  track  in  Figure  11.  Near  range  compression  due  to  the  slant  range 
presentation  partially  masks  the  difference  in  layover  between  the  two 
iniages.  Unfortunately,  the  scan  lines  in  the  IDECS  display  interfere  with 
the  study  of  the  two  simulations,  so  they  must  be  viewed  at  a distance. 
However,  differences  can  still  be  distinguished  in  the  category  discrim- 
ination ability,  fading,  object  brightening,  shadow  and  layover.  The  center 
pyramid  in  the  Figure  jg  has  many  cells  located  in  bins  along  a line 

parallel  to  the  flight  track,  whereas  Figure  11  has  a larger  area  of 
brightening  (lower  intensity)  constituting  a triangular  section. 

In  conclusion,  it  has  been  shown  that  the  SEAR  simulations 
correctly  represent  layover  and  shadow.  The  use  of  the  artificial  data 
base  makes  possible  the  quantification  of  these  radar  phenomena  for  vali- 
dation of  the  SCAR  model.  Empirical  backscatter  data  at  the  frequencies 
and  polarizations  of  Figure  7 lends  credence  to  the  idea  that  the  SLAR 
simulation  is  very  valuable  for  system  optimization  tests,  especially 

IDECS  - Acronym  for  jjnage  D^i  sc  r im  i na  t i on  , Enhancement,  Combination  and 
Rampling,  (an  image  processing  stationT,  Remote  Sensing  Labora- 
tory . 
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in  the  area  of  backseat  ter  category  discrimination.  It  is  not  possible 
to  make  statements  about  the  overall  adequacy  of  the  closed  system  im'>de  1 
for  it  is  not  yet  possible  to  objectively  measure  the  quality  of  the  simu- 
lated images  and  their  suitability  for  specific  applications.  Future 
studies  are  being  aimed  toward  the  goal  of  defining  image  quality  fac- 
tors for  very  specific  applications,  but  until  these  factors  can  be 
cjuantified,  the  value  and  quality  of  the  simulation  is  left  to  the  subjec- 
tive opinion  of  engineers  and  radar  interpreters.  It  should  be  noted, 

however,  that  simulations  of  real  scenes  produced  at  RSL  have  consistently 
been  mistaken  for  real  imagery  by  staff  radar  interpreters. 
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Cone  1 us  i f>ns  ant)  RecomnioncJat  ions 


The  ca  leu  I ah  If'  response  of  an  arlilicial  data  base  has  been  cjsed 
to  explore  the  accuracy  of  the  SLAR  simulation  programs  which  embody 
the  theory  briefly  described  in  this  document.  The  closed-system  model 
v-jhich  has  been  developed  to  imitate  ground-radar  interactions  has  been 
shown  to  have  considerable  flexibility  by  the  application  results 
presented.  The  simulated  images  shown  herein  illustrate  effects  such 
as  (1)  layover;  (2)  shadow;  (3)  frequency/ polar! zat ion ; (k)  near  range 
compression;  and  (5)  category  discrimination.  The  value  of  image 

simulation  for  pre-evaluation  of  future  radar  systems  has  been  stressed 
because  recent  developments  since  the  publishing  of  Reference  [2]  have 
shown  the  utility  of  a geometrical  data  base  (with  empirical  backscatter 
characteristics  and  computer  generated  shapes)  in  conjunction  with  the 
SLAR  simulation  programs  of  the  Remote  Sensing  Laboratory. 

The  lack  of  a large  catalogue  of  empirical  backscatter  data  (VV,  VH , 
HH,  HV,  circular  direct  and  cross  polarizations)  and  computer  resources, 
have  hampered  the  simulation  efforts  to  some  degree.  The  apparent 
advantages  of  accurately  predicting  system  performance  prior  to  construc- 
tion as  well  as  use  in  guidance  systems  would  seem  to  be  incentive  to 
apply  what  has  been  learned  about  radar-ground  interactions  through  use 
of  this  particular  closed-system  model. 

The  radar  image  simulation  methods  developed  will  next  be  applied 
to  a real  data  base  (Pickwick,  test  site)  to  measure  its  performance 
and  to  validate  the  model  by  comparison  of  simulations  with  real  imagery 
of  the  same  site.  Subsequently,  experiments  will  he  performed  at  the 
Engineering  Topographic  Laboratories  to  determine  the  suitability  of 
the  simulated  imagery  for  missile  navigation  systems  in  the  develop- 
mental stage. 


1 

( 

i 


Holtzman,  J.C..  V.H.  Kaupp,  R.L.  Martin,  E.E.  Komp,  and  V.S.  Frost, 
"Radar  Image  Simulation  Project:  Development  of  a General  Simula- 
tion Model  and  an  Interaction  Simulation  Model,  and  Sample  Results," 
TR  23^-13,  Remote  Sensing  Laboratory,  The  Utiiversity  of  Kansas, 
February  1976. 


TR  319-8 


APPENDIX  I 

The  SLAR  simulation  program  is  contained  in  this  section.  The 
data  which  must  be  supplied  to  the  program  were  previously  listed  in 
Table  I.  Due  to  interdependence  of  some  of  the  Radar  and  Ground 
matrix  values,  the  following  formulas  are  useful; 

••■Gr(9)  = R(8)  • Y 

*Gr(l0)  = R(6)  • Y • R(9) 

*Gr(ll)  = Gr(l0)/Gr(2)  + .5  1 1 Gr(ll)  ^ 5 

*Gr(l2)  = G(5)/Gr(ll) 

*Gr(13)  = R(7)/R(8) 

*R(15)  = R(16)  • Y 

The  asterisked  quantities  are  computed  within  the  body  of  the  program. 
When  changing  R(6),  for  example.  R(9)  must  be  altered  to  compensate 
such  that  Gr(lO)  has  the  desired  value. 

Parameters  for  Figures  4-5  and  4-6  are  given  below  for  a more  com- 
plete description  of  the  radar  systems  being  simulated  by  the  SLAR  pro- 
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CONSTRUCTION  OF  A GEOMETRIC  DATA  BASE 
FOR  RADAR  IMAGE  SIMULATION  STUDIES 


The  followinq  technical  report  (TR  3'9*1), 
prepared  by  the  Center  for  Research,  Inc., 
University  of  Kansas,  is  included  in  this 
report  to  provide  the  rationale  for  the 
discussion  of  Section  2. 
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ABSTRACT 


A computer  program  to  produce  a digital  ground  truth  data  base  of  a 
controlled  scene  of  geometric  solids  has  been  developed  for  radar  image 
simulation  studies  The  data  bases  which  can  be  made  by  this  program 
have  well-defined  geometric  properties  and  have  well-defined  boundaries 
separating  different  microwave  reflectivity  categories.  Since  the  geo- 
metric properties  of  the  features  in  a scene  are  of  paramount  importance 
to  radar,  a test  data  base  with  known  geometry  is  crucial  in  radar  image 
simulation  work.  Also  important  to  radar  image  simulation  work  is  the 
capability  to  easily  change  the  reflectivity  categories  of  the  various 
regions  contained  in  a scene.  All  of  these  features  are  available  in 
the  data  bases  which  can  be  produced  by  this  program.  These  data  bases 
are  ideally  suited,  therefore,  to  radar  image  simulation  studies;  they 
have  a known  input  to  the  simulation  program  and  produce  a calculable 
response. 


' ms  .■xr.^t-T-aaCa#»Sr~‘' 


INTRODUCTION 

A software  package  to  produce  a digital  ground  truth  data  base  of 
a controlled  scene  of  geometric  shapes  and  solids  has  been  developed  for 
radar  image  simulation  studies.  The  data  base  produced  by  this  computer 
program  provides  to  the  sinujlation  program  a known  input  with  well-defined 
interactions  and  relationships  between  the  various  features  contained  in 
the  scene.  The  uncertainties  in  the  geometry  and  dielectric  properties 
of  the  features  in  a scene  inherent  in  digital  data  bases  of  real  terrain 
are  removed.  Thus  the  data  bases  which  can  be  developed  from  this  program 
are  ideally  suited  to  provide  a known  calibrated  input  to  radar  image  simu- 
lation studies  of  an  investigative  nature;  they  have  a calculable  response, 
and  hence  can  be  thought  of  as  a standard  scene- 

The  most  difficult  aspects  of  the  radar  image  simulation  problem  are 
to  obtain  the  geometry  of  the  ground  features  and  knowledge  of  the  micro- 
wave  reflectivity  properties  of  the  ground  scene.  Some  of  the  more  common 
sources  of  these  data  are  high  resolution  aerial  photography,  maps,  and 
other  radar  imagery  at,  perhaps,  a different  wavelength  and  polarization. 
Independent  of  the  ground  scene  data  source,  it  is  necessary  to  extract 
the  desired  information  and  convert  these  raw  data  into  a suitable  digital 
data  base  (ground  truth  data  base)  for  input  into  a digital  computer. 
Typically,  the  scene  for  which  it  is  desired  to  simulate  the  imagery  of  a 
specific  radar  will  reguire  the  generation  of  a large  number  of  points  in 
the  ground  truth  data  base  For  example,  if  the  radar  system  to  be  modelled 
for  the  production  of  simulated  imagery  is  only  capable  of  resolving  fea- 
tures (objects)  which  are  separated  by  100  feet  or  more,  a 10  mile  by  10 
mile  scene  would  contain  approximately  279,000  independent  elements.  (50 
foot  resolution  would  reguire  more  than  a million  elements.)  Production 
of  simulated  images  of  scenes  of  a reasonable  size  can  require  the  develop- 
ment of  a very  large  ground  truth  data  base. 

In  a broad  sense,  there  are  two  classes  of  radar  image  simulations 
that  can  be  produced;  the  first  is  the  class  of  simulations  where  the 


simulated  radar  iinai|«‘ry  is  the  ((oa  1 , and  the  second  is  the  class  where 
the  imagery  is  not  the  goal.  The  simulations  of  the  first  class  are  pro- 
duced, typically,  to  reproduce  faithfully  the  radar  imagery  character- 
istics of  a particular  scene  based  on  the  real  properties  of  that  scene. 

For  this  class  of  simulations  it  is  necessary  to  generate  the  appropriate 
ground  truth  data  base  for  each  scene.  The  applications  of  the  simulated 
imagery  produced  in  the  second  class  are  unrelated  to  the  scene  itself. 

Some  examples  of  simulation  where  the  radar  image  is  not  the  goal  are 
found  in  troubleshooting  software,  in  sensitivity  studies,  in  parameter 
optimization  studies,  etc.  A standard  calibrated  data  base  with  known 
input  characteristics  would  be  ideal  for  this  class  of  simulation. 
Typically,  however,  an  existing  data  base  produced  for  some  other  appli- 
cation is  used.  This  approach  has  the  problem  that  it  is  difficult  (if 
not  impossible)  to  determine  the  origin  of  the  results  produced  in  the 
simulation  (e.g.  anomalies  might  be  caused  either  by  the  data  base  or  the 
study).  In  fact,  experience  has  shown  that  using  existing  data  bases  to 
troubleshoot  new  software,  to  determine  thresholds  of  sensitivity  for 
various  simulation  parameters,  to  evaluate  discrimination  capabilities, 
or  for  other  similar  studies  is  inefficient,  can  mask  results,  and  can 
lead  to  incorrect  interpretations  of  results.  A controlled  data  base  of 
known  geometry,  dielectric  categories  and  features,  and  having  a calculable 
response  is  reguired  for  studies  of  this  kind. 

Several  geometric  data  bases  have  been  produced  by  the  program  des- 
cribed here.  A specific  example  is  shown  in  Figures  1,  2,  and  3.  Figure 
I shows  an  isometric  line  drawing  of  this  data  base.  This  figure  illus- 
trates the  relative  orientation  of  the  various  geometric  shapes  included. 
Figure  2 is  a plan  view  of  the  same  data  base  showing  the  relative  heights 
and  positions  of  the  various  solids.  Figure  3 is  also  a plan  view  of  this 
data  base,  but  this  figure  illustrates  the  microwave  reflectivity  category 
assignments  made  for  this  one  data  base.  This  data  base  has  been  used  to 
troubleshoot  the  image  simulation  software  and  to  verify  that  layover, 
shadow,  range  compression,  local  angle  of  incidence,  and  fading  were  all 
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Figure  3-  Microwave  Reflectivity  Category  Assignments 


properly  simulated.  It  has  also  been 
effects  of  chanqinq  carrier  frequency 
scene.  It  is  expected  that  this  data 
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base  will  be  invaluable  in  future 


METHODOLOGY 


The  software  design  for  the  construction  of  the  various  geometric 
solids  was  made  as  simple  as  possible.  The  various  geometric  solids  which 
can  be  selected  for  inclusion  in  a data  base  to  be  constructed  by  this  pro- 
gram are  available  as  subroutines.  The  objects  available  as  subroutines 
and  the  options  for  user  specification  are  listed  in  Table  1.  This  pro- 
gram creates  a grid  matrix  of  an  arbitrary  size.  (Present  operation 
creates  1000  x 1000.)  Cell  size  is  user-specified,  but  must  be  square. 

Within  this  grid  matrix,  3-dimensional  representations  of  regular 
surfaces  are  constructed  by  user  specifications  (details  of  available 
structures  to  follow).  In  addition,  a specific  microwave  reflectivity 
category  may  be  assigned  to  each  object.  The  elevation  and  category 
assignment  (i.e.  reflectivity  category)  for  each  cell  is  packed  into  one 
word;  the  rightmost  six  bits  are  reserved  for  category  number  (giving 
possible  range  of  values  0 ^63).  and  the  remaining  upper  bits  are  reserved 
for  the  elevation.  Those  cells  where  no  object  is  defined  are  implicitly 
defined  with  0 elevation  and  0 category. 

This  computer  program  will  produce  data  bases  as  large  (or  as  small) 
as  is  desired.  To  prevent  problems  caused  by  computer  core  limitations, 
the  data  base  is  constructed  in  horizontal  strips.  In  the  first  pass, 
all  the  columns  of  the  first  N rows  of  the  grid  matrix  are  constructed 
and  stored  on  an  external  device;  pass  2 produces  the  next  N rows;  pass 
3 the  next  N rows,  etc.,  until  the  matrix  is  complete.  N is  chosen  by 
the  user  so  that  (N  ROWS  * M COLUMNS)  is  a "reasonable"  amount  of  core 
requirements.  For  the  specific  data  base  herein  described  (Figures  I,  2, 


TABLE  I 

Geometric  Solids 


Solid  Opt i ons 


Cube  User  supplies  start  row  and  start  column,  length  and  width 

of  base,  the  angle  of  rotation  in  degrees,  and  a microwave 
reflectivity  category  number  assignment. 

Wedge  User  specifies  start  row  and  start  column,  length  of  wedge, 

maximum  height  of  the  wedge,  slope  of  the  front  face,  width 
of  the  flat  top,  slope  of  the  back  face,  angle  of  rotation, 
and  a microwave  reflectivity  category  number  assignment. 

Pyramid  Produces  a pyramid  of  arbitrary  height  with  a square  base. 

User  specifies  start  row  and  column,  length  of  base,  angle 
of  rotation,  and  microwave  reflectivity  category  number. 
(Slope  of  the  sides  is  implicitly  determined  by  the  height 
and  length  of  the  side.) 

Hemisphere  User  specifies  row  and  column  of  center,  radius,  and  micro- 
wave  reflectivity  category  number. 

Tower  Produces  a cylinder  standing  on  end.  User  specifies  start 

row  and  column  of  center,  radius,  length,  and  microwave 
reflectivity  category  number. 

Cylinder  Produces  half  of  a cylinder  lying  on  its  side.  User  speci- 

fies start  row  and  column  of  center  line,  radius,  length, 
and  microwave  reflectivity  category  number. 

Patch  Assigns  to  a given  number  of  matrix  cells  the  desired 

microwave  reflectivity  category  number  with  no  elevation. 

Checkerboard  Assigns  to  a series  of  matrix  cells  each  of  the  different 
microwave  reflectivity  categories  consecutively  with  no 
e 1 evat ion . 

Note: 

Patch  and  checkerboard  are  available  to  compare  changes  in  the  return 

strength  of  the  various  categories  at  different  incident  angles.  Patch 

may  provide  a reference  greytone  next  to  an  object  to  analyze  the  effect 

of  changing  the  local  incident  angle. 


nnd  3),  a 1000  x 1000  matrix  was  produced  in  20  passes.  Each  pass  pro- 
duced a section  of  50  rows  x 1000  points  which  recpiired  a single  block 
of  50K  words  of  core  (compared  to  lOOOK  which  would  be  required  if  the 
entire  matrix  were  constructed  in  one  pass).  The  multiple  pass  concept 
does  not  sacrifice  execution  speed  since  only  one  object  can  be  constructed 
at  a time.  And  since  tfie  outi)ut  is  performed  on  entire  lines,  the  data 
base  is  built  as  a single  unit,  eliminating  the  complexities  of  mosaicking 
sub-images  together.  The  one  sacrifice  that  is  made,  however,  is  that 
the  objects  produced  must  lie  completely  within  a single  strip;  an  object 
cannot  be  constructed  that  begins  someplace  within  one  horizontal  strip 
and  overlaps  to  some  point  within  another  strip.  Aside  from  this  minor 
limitation,  the  computer  program  is  the  most  general  in  design  it  was 
possible  to  produce. 


FLOW  CHART 


The  computer  program  (geometric  data  base  construction)  has  been 
designed  to  be  as  versatile  as  possible.  The  size  of  the  data  base  and 
the  size,  location,  orientation,  and  shape  of  the  various  geometric  solids 
are  easily  controlled.  Figure  ^4  presents  the  basic  operations,  in  flow 
chart  form,  of  this  computer  program.  Digital  data  bases  of  geometric 
solids  are  constructed  according  to  the  user's  design  and  specification. 

Also  shown  is  a listing  of  the  eight  geometric  shapes  available  as  sub- 
routines in  this  prop ram.  A copy  of  the  program  is  included  in  Appendix 
I of  this  work. 

To  use  this  program,  it  is  necessary  first  to  design  the  desired 
data  base.  The  size,  location,  orientation,  and  shape  of  each  geometric 
shape  must  be  specified.  The  computer  program,  as  listed  in  Appendix  I,  will 
produce  a digital  data  base  matrix  constructed  of  1000  elements  by  50 
rows  at  a time  for  as  many  repetitions  as  are  desired.  (A  matrix  of 
1000  X 1000  data  values  would  require  20  repe t i t i <uis . ) If  it  is  desired 
to  construct  a data  base  of  a different  size,  the  dimension  statements 
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must  be  changed  from  the  values  specified  (specified  in  the  program  listed 
in  Appendix  f ''  to  the  desired  values.  The  only  real  constraints  on  data 
base  size  are  imposed  by  the  core  size  of  the  computer  used  to  construct 
the  data  base.  After  the  data  base  size  has  been  defined  and  the  dimen- 
sion statements  of  the  program  have  been  specified  correctly,  the  data 
base  can  be  constructed.  Reference  to  the  flow  chart  shown  in  Figure  A 
will  make  clearer  the  following  discussion.  For  each  geometric  shape  to 
be  included  in  each  horizontal  strip,  two  data  cards  are  required  to  be 
input  to  the  program.  The  first  card  is  a number  between  0 and  9 that 
specifies  the  surface  to  be  constructed,  and  the  second  card  specifies 
the  required  construction  data.  (These  requirements  are  listed  both  in 
Table  1 and  in  the  various  subroutine  listings  in  Appendix  I.)  Two 
cards  are  required  for  the  program  for  each  shape  desired  in  a strip.  A 
data  card  of  0 signifies  that  no  more  geometric  solids  are  desired  for  a 
given  strip  and  that  construction  is  to  begin  for  the  next  strip.  In  this 
way,  the  data  base  is  constructed  strip  by  strip  until  the  data  base  is 
completed . 

The  data  bases  which  can  be  constructed  from  this  program  are  not 
likely  to  be  encountered  in  actual  radar  use.  But  they  can  be  very  use- 
ful for  testing  radar  simulation  programs  and  problems  since  all  aspects-- 
height,  local  slope,  category  ass i gnment--are  precisely  known.  Further- 
more, any  or  all  of  these  parameters  can  be  changed  by  the  user  in  a very 
controlled  manner  to  study  the  relative  effects  and  importance  of  each, 
as  well  as  to  evaluate  the  simulation  techniques.  These  changes  can  be 
accomplished  either  by  constructing  a new  data  base  for  each  arrangement 
or  by  incorporating  this  program  in  the  simulation  model  to  produce  the 
desired  effect.  One  can  also  study  the  importance  of  the  size  of  the 
data  base  cells  by  changing  the  appropriate  parameter  during  the  pro- 
duction of  a data  base.  There  are  many  more  potential  uses  for  a data 
base  like  this,  and  to  increase  the  value  of  this  program  by  making  the 
output  more  generally  applicable,  facilities  to  rotate  the  data  base  to 
run  radar  simulations  with  different  flight  paths  have  been  included. 
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Listing  of  Computer  Program  to  Construct 
Geometric  Data  Bases 
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reflectivity  for  input  to  a radar  image  simulation 
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C THI^  PKOGKAIV  CHtATFS  AM  ARTIFICIAL  n A I A RASE  f'TMSISTIMG 
C CF  REGULAR  S-'DIMEMSICMAL  GEOA’FTRICAI,  shapes  AS  SPECIFirf) 

C EY  THE  USER 
C 

C IHF  DATA  EASE  IS  CONSTRUCTEF)  IM  RFfTANELLAR  FOR'^'AT. 

C EACH  CELL  REPRESENTS  A FIXED  SOUARF  AREA,  THE  ICMGTH  OF 
C THE  SIDE  OF  EACH  CEll  IS  OFTFRMINFO  nv  THE  VARIARLE  — 

C CELS  W 

C EACH  CELL  IS  CHARACTERIZED  RY  A\  ELEVATION  A Ml)  A CATEGORY 
C THIS  INEORyATICM  IS  FAfKFr*  INTO  0’'=^  '.vCCD. 

C the  low  ORDER  6 4ITS  ARE  RESERVED  FOR  THE  CATEGORY^ 

C THE  HIGH  ORDER  HITS  FOR  TH'^  FIEVATIPN 
C 

C TO  USE  THIS  PROGRAM  THE  USER  SPECIFICS  THE  OGJFCTS  DESIRED 
C ALCNG  WITH  THE  NECESSARY  DESCRIPTIVE  PARAMETERS  AND  THE 
C CCCRCINATES  CF  THF  ICWER  LEFT  CCRNFR  FOR  PLACE'^FNT  IN  THF 

C CATA  BASE  ( OR  THE  CENTER  FOR  CYLINDRICAL  AND  SPHERICAL 

C OBJECTS)  PLUS  THE  DESIRED  CATEGORY 
C 

C ALL  UNSPECIFIED  CELLS  HAVE  IMPLICIT  ELEVATION  A D CATEGORY  o 
C 

C TO  MINIMIZE  CORE  REQUIREMENTS  THE  DATA  EASE  IS  CONSTRIJCTED 

C IN  STRIPS  50  CELLS  WIDE  (THIS  NUMBER  MAY  BE  VARIED  BY 

C ALTERING  the  DIMENSION  STMtT).  THE  FIRST  FIFTY  ROWS  ARE 
C CONSTRUCTED  AND  PUT  ON  TAPE  WHEN  CO^'OCFTF.  THEN  THF  NFXT 
C FIFTY  ROWS  ARE  CONSTRUCTED  AND  PUT  TAPE»  AND  SO  ON 

C THF  USER  SPECIFIES  THE  COMPLETION  OF  A STRIP  ON  INPUT 
C FY  IAPUTTING  0 FOR  CHOCSF 

C THEREFORE  AN  ARBITRARY  NUMBER  CF  ROUS  MAY  BE 
C CONSTRUCTED.  THE  NUMBER  OF  COLUMNS  IS  LIMITED  BY  THE 
C dimensions  of  the  array  MATRIX  (IN  THIS  CASE  IPPP) 

c 

c CJECTS  may  not  cross  OVER'  THE  LVO  UNDARIFS  OF  STRIPS 
C EACH  OBJECT  specified  «'(JST  FIT  CC^'El.  ETFl.  Y INSIDS  the 
C STRIP 
C 
C 

implicit  integer  (a-W) 

c 

c the  following  variables  are  REQIIIRFd  py  most  or  all  of 
C The  subroutines  and  so  are  placed  in  blank  common  and  so 

C DC  NCT  HAVE  TO  HE  PASSED  AS  FORA'AL  PARAMETERS 
C 

C E>UF  - small  I/O  buffer  FOR  DEBUG  WRITE  STMT  IN  VARIOUS 
C SUBROUTINES 

C MATRIX  - THAT  PORTION  OF  THE  DATA  pasE  PPiNr-  CONSTRIJCTED 
C THAT  IS  PHYSICALLY  IN  CORE 

C CLLSIZ  - SIZE  CF  EACH  CELL.  IN  DATA  daSE  (IN  FEET) 

C STROW^STCOL  - STARTING  ROW  AND  C Cl  U^'N  OF  THF  SUPFACF  TO  BE 


I 
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C CONSTRUCTED  (FOR  PFCTAMGULAR  CRJFCTS  THIS  IS  LOWER  LEFT 

c corner;  for  cylindrical  objects  the  center) 

C ***N0TE***  STRCW  NUST  be  given  f'COULC  50  (SINCE  ONLY  A 

C 5C  ROW  STRIP  OF  THE  ENTIRE  DATA  PASF  IS  IN  CORE 

C AT  ANY  ONE  T I E > 

C CAT  - CATEGORY  ASSIGNMENT  FOR  THE  SURFACE 

C ANC-  = ANGLE  OF  ROTATION  (APPLICABLE  ONLY  TO  RECTANGULAR  SHAPES) 
C THE  OBJECT  WILL  RE  ROTATED  COUNTER-CLOCKWISE  ABOUT  ITS 

C LOWER  LEFT  CORNER  THE  SPECIFIED  NUMBER  nF  DEGREES 

C 
C 

COMMON  8UF(20>/  matri X ( lono/ 50 /CElsi z »stcol/S tpow/C a T/ an f 
DIMENSION  LINE(500) 

c 

CELSIZ  =20 
HALFCEL  * CELSIZ/2 
C 

00  620  1*1/70 
DO  62C  JS1/50C 
62C  MATRIX{J/I)=10 
C 
C 

C CONSTRUCTION  OF  EACH  SHAPE  REOUIRES  TWO  DATA  CARDS 
C FIRST  CARO  IS  A NUMBER  BETWEEN  0 AND  9 (I?  FOR'^AT) 

C THAT  SPECIFIES  THE  SURFACE  TO  BE  CONSTRUCTED 

C C - BEGIN  NEW  STRIP  OF  DATA  BASE 

C 1 - PATCH 

C 2 - CUBE 

C 3 - WEDGE 

C 4 - FRYAMIO 

C 5 - HEMISPHERE 

C 6 - CYLINDER 

C 7 - TOWER 

C 8 8 - CHECKERBOARD 

C 

C SECOND  CARO  NECESSARY  DESCRIPTIVE  PARAMETERS  FOR  THAT 
C OBJECT  --  SEE  SPECIFIC  ROUTINE  FOR  DETAILS 
C 

C *«*NCTE**»  STROW  and  STCOL  ARE  TC  BE  SPECIFIES  IN  TERMS 
C CF  MATRIX  CELLS  (ROW/COL)  - ROW  MODULO  50  - 

C FOR  EASY  PLACEMENT  OF  OBJECTS  RELATIVE  TO  ONF  ANOTHER 
C ALL  OTHER  PARAMETERS  ( L E NG T H / W 1 D T H / HEIGHT)  ARE  TO 

C BE  GIVEN  IN  FEET  SO  THE  USFR  CAN  EASILY  OESCRIFIF  THF 
C PHYSICAL  SIZE  DESIRED 
C 
C 

1 READ(O5/C5/EN0=9C0)  CHOOSE 

5 fORMAT(l2) 

C 

IFICHOOSE  .EQ.  C)  GO  TO  SGf 

GO  TO  (10»20/3r#40»5n/6n,7o,?C)/  CHCnSF 


C PAT  INI' I I (UTFUT  FFt<oi.'  F S A ( ( Arr  D F A p f * r n 

C 

UR  1 r E (6» >»  1 D chccsr 

81  C fORRAT(//#'  »**UARMA4C.  tAllwnprt;  VA(.  (IF  TO 

1 SPECIFY  NEW  OPJECT.  IT  ('A<;  f- F F ►,  IGNOErn'//) 

C 

C SKIP  CARC  WITH  PARAVETFRS  FOR  PAD  Oi.-JEFT 
C 

RtAn(C5#C5/FNO  = <irC)  CHOOSf 
CO  TO  1 

c 

c 

C PATCH 

c 

C PATCH  - SQUARE  AREA  OF  OATA  PACr  a c c;  j r,  A F 0 FPECIFlPr)  caTFOORY 
C NO  ELEVATION  IS  ASSIGNEO 

C 

C SIZE  - DESIRED  LENGTH  CF  SIDE  FOR  Thf  SOUARF  ARFA 
C 

c 

1C  READ(05/11)  STROw,STCOL #ST;r,r at 
11  FOHI^ATdOIS)  ’ 

CALL  PATCH  (SIZE) 

CONTINUE 
CO  TO  1 
C 

C C U P F 

C 

C CUBE  - rectangular  SHAPE  OF  ARPTRAy  F 0 N S T a EirvATlOM 

c 

C LEN  - length  in  X OIPFCTION 
C WIC  - WIDTH  IN  Y DIRECTION  (VtiTTKAL) 

C HGT  - HEIGHT  CONSTANT  FLFVATION  T^  AFSICafO  rpLL 

C 

C 

2C  KEAD(05#11)  S TRCU/ S T col / ANC, L E N# W I 0 /HG T /C a t 
CALL  CUBE ( Lf N, U I D/ HGT ) 

CONT INUE 
GO  TO  1 
C 

C WEDGE 

C 

C WEDGE  - OBJECT  WITH  SLOPING  FRONT  fafF/  a f|ST  TOP 

C OF  CONSTANT  ELEVATION  AND  ARPITRARY  WIDTH  AND  A SLOFIa.’O  Fach 

C FACE 

C 

C LEN  - LENGTH  OF  WEDGE  IN  X OIRFCTIPT  (AfROSS) 

C SLOPE  - SLOPE  CF  FRONT  FACF 
C HGT  - MAXIMUM  FLFVATION  OF  WFDGF 

c TCP  - width  of  flat  top 

C PSLOFE  - SLOPE  OF  BACK  FACE 
C 


3C  REA0(C5/11)  STROW/STCOL^ANG/LEN/SlCFE/HGT^TOP/OSLOPE# CAT 
call  wedge (LEM/SLO PE »HGT»TOP»eSLOPF) 

CONT  INUE 
CC  TO  1 
C 

C PY  R AM  D 


C 

C FRYAMO  - REGULAR  FOUR  SIDED  PYRAMID 
C 

C LE^  - SIZE  OF  BASE  OF  PYRAMID  (WILL  HAVF  SOUARP  PA<;f) 

c HG T - height  of  PYR am D 

C 

C 

4C  READ(05#11)  STROU/STCOL/ANG/LC'/HGT/CAT 
call  PRYAMID (LEN^HGT) 

CONTINUE 
GO  TO  1 
C 

C HEMSPHEPE 

C 

C HEMSPHERE  - HEMISPHERE  WITH  CENTER  AT  ( S T ROW  / S T C OL  ) 

C 

C RAC  - RADIUS 

C 

C 

5C  REAO(C5/n)STROW#STCOL/RAO#CAT 
CALL  HE«ISF(RAD) 

CONTINUE 
GO  TO  T 
C 

C CYLINDER 

C 

C CYLINDER  - CYLINDER  LYING  ON  ITS  SIDE 
C 

C RAD  - RADIUS  OF  CYLINDER 

C LEN  - LENGTH  OF  CYLINDER  (IN  X DIRECTION) 

C 

C 

6C  fiEAD(05#11)  STR0W,STC0L/RAD/LEN/C AT 
CALL  CYL(RAD,LEN) 

CONT INUE 
GO  TO  1 
C 

C TOWER 

C 

C tower  - cylinder  standing  ON  END 

c 

C RAO  - RADIUS  OF  CYLINDER 
C HGT  - HEIGHT  OF  CYLINDER 
C 
C 


I 


c 


/(I  AO  (05/  1 1 ) STRCW/STCOL  /RAO/HGT/f  AT 

gall  T0WER(R AO/HGT ) 

CONTINUE 
CO  TC  1 
C 

C CHEC  KER"CARD 


c 

C CEFCKtRHOARO  - ROW  AND  COLO^’N  Of  f- A I r H F <;  P F P r F <;  f N T I N’ f,  THF 
C VARIOUS  CATEGORIES  (FOR  CONCARISCN  OF  REi-vyiVE 

C GRFYTCNES  FOR  THF  CATFGORIFS) 

C 

C 

C 

fiC  READ(G5/11)  STRCW/STCOL 
CALL  CHKRBD 
CONTINUE 
CO  TC  1 
C 

c 

C OUTPUT  completed  STRIP  OF  DATA  PAFF  TO  tapf 
C IN  BINARY  UNFORMATTED  FORM 
C 

5fC  CO  6C0  Ia1/50 

WRITE(Cl)  (MATR I X( J, I ) / J= 1 , 1 CrO) 

6rc  CONTINUE 
C 

C RESET  ALL  MATRIX  CELLS  TO  0 HFFORF  UPGINNING  NFxr 
C STRIP  OF  DATA  EASE 
C 

DO  61C  1=1/50 
DO  610  J=1/10C0 
610  MATRIX(J/I)  = 0 
GO  TO  1 
9rC  STOP 
END 

C SUBROUTINE  PATCHCSWE) 

SUBROUTINE  PATCH (SIZE) 

C 

IMPLICIT  INTEGER  (A-W) 

COMMON  BUF(20)/  A R R A Y ( S PC’ / 70  ) / C F L S I Z / S T C OL  / S T R /C  A T / ANG 
C 

c 

C SIZE  DIMENSION  CONVERTED  TO  NUMBER  OF  •'ATRIX  CELLS 
C 

SIZE  = SIZE/CELSIZ 
C 

ro  ICC  1 1 = 1/ s iz E 
I s I 1-1 

CO  ICO  JJ»  1/ S I Z E 
J sJJ-1 

c 


1 


-JBW-XS. 


C /IFFRCPRIATE  CELLS  ARE  ASSlGN'Fn  THE  C F <;  I G N A T E 0 CATEGORY 
C 

100  ARRA V( STCOL* I »S TPOW+J ) = CAT 
RETURN 
END 

C SUBROUTINE  C UB E ( L E N# W I 0 / H G T ) 

SUBROUTINE  C UB E ( L E N / W I D / HG T ) 

C 

IMPLICIT  INTEGER  (A-W) 

COMMON  BUF(20)»  ARRAY(5nO#7n)/CELSIZ,STCOL»STROW/CAT/ANG 
DATA  LC#BC»LR/PR/500/0, 500/0/ 

C 

C 

I = -1 
C 

HALFCEL  = CELSrZ/2 
C 

C ANGLE  OF  ROTATION  OF  BASE  OF  CUBE 
C 

ZANG  = FLOAT (AN6)/57. 295 
C SINE  OF  ANGLE  CF  POTATION 
ZSIN  s SINCZANG) 

C COSINE  OF  ANGLE  OF  ROTATION 
ZCOS  = COS(ZANG) 

C 

C 

C 

C THESE  DC  LOOPS  CALCULATE  THE  CELLS  THAT  WOULD  BF  INCLUDED 
C IN  THE  CUBE  IF  THERE  WERE  NO  ROTATION 
C IF  THE  OBJECT  IS  TO  BE  ROTATED  (ANG  =/  C) 

C A ROTATION  OF  AXES  IS  EFFECTED  BY  THE  TRANSFORMATION 
C XNEW  = X A COS(ANG)  - Y * SINCANG) 

C YNEU  » X * SIN(ANG)  ♦ Y * COS(ANG) 

C 

C THEN  THE  CELLS  REPRESENTED  BY  XNFW  AND  YNFW  ARE  ASSIGNED 
C the  APPROPRIATE  VALUES 
C NOTE  IF  ANG  = C NO  CHANGE  IS  MADE 
C 

C THIS  SAME  technique  IS  USED  FOR  ROTATION  OF  ALL  OBJECTS 
C 

CO  ICO  IV  = HALFCEL/WID/CELSIZ 
Y = lY 

YSIN  = ZSIN  A Y 
YCOS  * ZCOS  A Y 

DO  ICC  IX  « HAL FCEL/LEN/CELS  I 7 
X » IX 

COL  * ZCOS  A X - YSIN 
ROW  « ZSIN  AX*  YCOS 
COL  * COL/CELSIZ  ♦ STCOL 
ROW  « ROW/CELSIZ  ♦ STROW 
C 

C LR/  LC  - LEAST  ROW  AND  LEAST  COLUMN  INCLUDED  IN  THIS  OBJECT 


I 

j 


C en,  0C  - 3IGGEST  row  AND  COLU'^N  I N r 1. 1)  (i  (=  D 

C THESE  PARAMETERS  REQllIPro  fCP  The  SUPPOUTINF  Fill  IN 

C 

IFICCL  .LT.  Lf)  LC  = COL 
IF  (CCL  .GT  . HC)  BC  = COL 
IFCEOW  .LT.  LR)  LR  = POW 
I F ( ROW  .GT.  BP)  PR  = POW 
C 

C OEPUGWRITFSTATEMFNT 

C 

GC  TO  90 

1=1+2 

BUF  ( I )=  COL 

BUF  ( 1 + 1 ) = ROW 

IF  { 1 . L T . IS  ) GC  TO  90 

URITe(6/60)  (PUF(M), 

60  FORM  AT (8 (5 X, 21 S ) ) 

I=-1 

90  f L 0 ( 18, 1 2/ ARRA Y ( COL/ROW ) ) = F L D ( ? A , 1 ? , HG T ) 

C 

C PROPER  ELEVATION  AND  CATEGORY  IS  ASSIGNED  TO  EACH  OF  THFSp 
C CELLS 
C 

ICC  FLC(  3C,6,ARR AY( COL,ROW) ) = F L D ( 3 0 , 6 , C A T ) 

C URITE(6/6r))  (PUF  (M)  ,M=1  , 1 6) 

DO  1C  1 = 1,16 
1C  EUE(I)=n 

CATEG  = CAT 


C SINCE  ROTATION  WAS  DCNF  ON  DISCRETE  POINTS  IT  IS  POSSIBLE 

c THAT  There  are  some  holes  (cells  into  which  no  potateo 

C POINT  WAS  PLACED)  IN  THIS  OBJECT 
C 

C SUBROUTINE  FILLIN  CHECKS  OBJECTS  FOR  SLCH  'HOLES* 

C AND  FILLS  them  IN  WITH  THE  APPROPRIATE  ELEVATION  AND 
C CATEGORY 
C 

CALL  E ILLIN(LC,BC,LR,8R, CATEG) 

CONT  INUE 

RETURN 

END 

C Suo  routine  WEDGE(LFN,SLCPr,HGI,TO'’,FiSLOFF) 

SUBROUTINE  WEDGE  (LEN, SLOPE, HGT, TCP, ESLOPE) 

C 

implicit  INTEGER  (A-W) 

COMMON  HUE(20),  ARRAY(5nn,70),C'"LSIZ,STCOL,STROW,CAT,AKT, 

DATA  LC,BC,LR,eR/scn,o,src,n/ 

C 

kRITE(6,6S)  STCCL,STROW 
6S  FORMAT!/, 3 CX, 'WEDGE', Sx, 216,//) 


FT  =-2 

FALFCEL  = CELSI7/? 

C 

C angle  OF  ROTATION  - SFF  ClI'^F  FOP  OFTAIL?  OF  ROTATION 
C 

ZANG  = FLOATIAMrO/S/.^PS 
ZbIN  = SIN(ZANG) 
zees  = CCSTZANG) 

C 

ZZTOP  = FLOAT ( SLCPF ) /b 7 . P9S 
C 

C TANGENT  OF  SOP  FOR  FRON  FAC*-'  OF  WEPGF 
ZTAN  = SIN(ZZTOP)/COS(ZZTOP) 

CNT  =C 
C 

C ELEVATION  INITIALIZED  TO  ZERO 
ELEV  =0 
C 

C UICTP  OF  TOP  IN  NUMBER  OF  CELLS 
TOP  = TOF/CELSIZ 
INCR  =0 

lY  = HALFCEL-CELSI7 
ICC  ELEV  = ELEV+INCR 

IFCELEV  .LT.O)  GO  TO  200 
lY  = lY^  CELSIZ 
Y = I Y 

YSIN  = ZSIN  * Y 
YCOS  = ZCOS  * Y 

c 

CO  11C  IX  =HALFCEL»LEN,CELSI  Z 
X = I X 

COL  = ZCOS  * X - YSIN 
ROW  = ZSIN  * X + YCOS 
COL  = CCL/CELSI Z ♦ STCOL 
ROW  = ROW/CELSIZ  ♦ STROW 


PARAMETERS 

FOR 

F I LL  I N 

AS 

DESCPIPEO  IN  CURE 

IF (COL 

.lt 

. LC) 

LC 

s 

COL 

I F (COL 

.6T 

. SC  ) 

PC 

= 

COL 

I F ( ROW 

.LT 

. LR) 

LR 

s 

ROW 

IF (ROW 

• GT 

. BR) 

PR 

s 

ROW 

GO  TO  90 

FT  =PT-f5 

OEPUG  PRINT 

STATEMENT 

eUF (PT ) 

= COL 

euF (PT41 ) aROW 

EUF (PT*2)* 

ELEV 

I F (PT  . 

LT. 

13) 

GO 

TO 

qr. 

URIIt(6#60)  (RUF('')#  f^=1,1S) 

6C  FOKVAT(5(5x,?IS,?X,IS)) 

FT  =-? 

F L D ( 1 8# 1 AR R A Y ( COL»ROU ) ) = F L C ( ? A , 1 ? , F L F V ) 

FLD( 30»6#ARRAY{ COL/ROW)  ) = F L D ( ^ r , Y , C A T ) 

11C  CONTINUE 


C CHECK  FOR  REACF-INC,  TCF^  OF  WFOr.  F 
C 

IF  (ELFV  .GE.  H'-T)  GO  TO  1 0 0 
C 

C HAVE  NOT  REACHED  TOP  SO  INCREMENT  FIEVATION  ^00  NEXT  ROW 
C PFiCCEEO  TO  NEXT  ROW 


INCR  = Z TAN  * FLCA  TCCEL  SI  7 ) 

GO  TO  TOO 
C 

150  CN7  = CNT  +1 

C HAS  flat  top  of  WEDGE  TEEN  COMPLETED'’ 

C 

IFCCNT  .GT.  TCP)  GO  TO  160 


C NO/  SO  NO  CHANGE  IN  ELFVATION  FOR  NEXT  ROW 
C 

INCR  sO 
GO  TO  100 

16C  IFCEjSLOPE  .GE.9C)  GO  TO  POr. 

73  = FLOAT  (BSLOPE)  7 5 7 . ?95 


C BEGIN  descending  ALONG  HACK  FACE  OF  WEDGE 
C MODIFY  SLOPE  FACTOR 

C NOTE  IT  IS  NEGATIVE  SO  THAT  ELEVATION  WILL  PE 

c decremented 
c 

ZTAN  = -SIN( Z3) /COS(23) 

INCR  = ZTAN  * CELSI7 
GO  TO  ICC 

2rc  continue 

C2CC  URITE<6/60)  { BU F ( M ) /M= 1 / 1 5 ) 

DO  11  1*1/15 
11  FUF  ( i)=n 

CATEG  * CAT 

CALL  FILLIN  (LC/RC/LR/BR/CATEG) 

CONT  INUE 
RETURN 

END  , 

C SUPKCUTINE  PRYAMID(LEN/HGT) 

SUBROUTINE  P R Y A M I D ( L E N / HG T ) 

C 

IMPLICIT  INTEGER  (A-W) 


COMMON  BUf(ZO),  /'RRA  y ( S f-n,  ) , OPLS  I ? / S r COL  TROW /C  A T , ANC, 
DATA  LC#BC,LR,HR/Srr,C,Sr.C,r/ 

c 

C ancle  Of  ROTATION  If  ANY  (SFF  CUBE  FOR  DETAIL!?) 

C 

ZANG  = FLOAT (ANG)/S7.?qs 
ZSIN  = SIN(ZANG) 

ZCOS  = C0S(ZAN6) 

C 

FALFCEL  = CELSIZ/^ 

C 

C HALF  THE  length  OF  THE  FINISHED  FRYAMIC 
C 

HFLEN  = LEN/2  + HALFCEL 
C 

C TANGENT  OF  SLOE  OF  PRYAMID  FACE 
C 

ZTAN  = fL0AT(H6T>/FL0AT (HELEN) 

C 

C 

C DO  LOOP  FOR  ROUS  OF  PRYAMID 

C BECAUSE  OF  SYMMETRY  FOR  EACH  POINT  RELCU  THE  CFNTERLINF 
C THERE  IS  A CORRE SSPCNO ING  POINT  EQUAL  DISTANCE  arOVE  IT 
C Y - IS  Y VALLE  FOR  LCUIER  CELL 
C Y2  - CCRRESSPONOING  CELL  ABOVE 
C 

CO  100  I YsHAL FCEL#HFLEN,CE LS  I Z 
Y = I Y 

YS IN  * ZSIN  * Y 
YCOS  = ZCOS  * Y 
Y2  = LEN  - I Y 
Y2SIN  = ZSIN  * Y2 
Y2CCS  = zees  * Y2 
C 
C 

C X VALUE  FROM  LEFT  TO  RIGHT 
L 

DO  150  IX  = HALFCEL^LEN^CELSIZ 
X » I X 

COL  = ZCOS  * X - YSIN 
ROW  * ZSIN  * X ♦ YCOS 
COL  = CCL/CELSIZ  hSTCOL 
ROW  = ROW/CELSIZ  ♦ STROW 
C 

C0L2  = ZCOS  * X - Y2SIN 
R0W2  = ZSIN  * X ♦ Y2C0S 
C0L2  = C0L2/CELSIZ  ♦STCOL 
fiOU?=  R0U2/CELSIZ  + STPOW 
C 

C parameters  for  F ILL  in 
C 


IF(CCL  .LT.LC)  LCsCOL 
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REMOTE  —ETC  F/6  17/9 

DAAK02-73-C-0106 

NL 


If  (CCl?  .1-1.  LC  ) LC  = COl  ? 

1 f (COL  . G f . HC ) 6C  = COL 
I F ( COL  2 .GT.  HC ) 8C  = COL  ^ 

1F(N0W  .LT.  LR)  LR  = ROW 
IF(R0W?  .LT.  LR)  LR  = ROW? 

1F(R0U  .GT.  BR)  ER  = ROW 
I F (ROW?  .GT.  BR ) HR  =ROW? 

C 

C IF  X < Y TH£^  elevation  AT  THAT  C F I,  L = XOI?»SLOPr  OF  FACE 
C (REPRESENTS  LEFT  FACE  OF  PYRAMID) 

F L E V = Z T A N * X 
IF  ( I Y . 6E  . IX)  GO  TO  7 P 
C 

C IF  X>Y  AND  Y<  LENGTH-XDIR  THEN  FLFVATICN  = YDIS*SL0PE 
C (FROM  FACE) 

C 

ELEV  = ZTAN  * Y 

IF(1Y  .GE.  (LEN  + CELSIZ-IX))  F.LFV  = ZTAN*(LFN'  + rFLSlZ-IX) 

C 

C ASSIGN  THE  TWO  SYMMETRICAL  POINTS  T K r APPRCPRIATF 
C ELEVATION  AS  JLST  CALCLLATFD 
C 

7C  fLO(18/12/ARRAY(COL#ROW))=  F L R ( ?4 , 1 ? , E L E V ) 

FLD(18#1?»ARRAY(C0L2/R0W?)  ) = Fi.  R(  ?4/1  ?,ELFi/) 

FLD(30^6»ARRAY(C0L,R0W))  = FLOISr^-E/CAT) 

FL  0 ( 3C»6/ ARRAY  ( CCL?#ROW?  ) ) = F L D ( 3 f # (S  , C A T ) 

15C  CONTINUE 
1CC  CONTINUE 

CATE6  = CAT 

CALL  FILLIN(LC/PC»LR/BR/CATFG) 

CONTINUE 

RETURN 

END 

C SUBROUTINE  H£MISF(RAO) 

SUBROUTINE  HEMISF(RAO) 

C 

IMPLICIT  INTEGER  (A-W) 

COMMON  BUF(20)/  array ( 5CC# 7n ) ,cflS I Z#S TC0L#R trow »C at, ANG 
F-ALFCEL  = CFLSI  Z/? 

PRAD  = RAO 
C 

C RADIUS  IN  NUMBER  OF  MATRIX  CELLS 
C 

NCELL  = FLCAT(PRAO)/FLOAT(CELSl/)  +.5 
/RAO?*  RAO**? 

C 

C CENTER  POINT 
C 

ARR A Y ( STCOL# S T ROW  ) = RAO 
C 

C UTILIZE  SYMMETRY  OF  HEMISPHERE 

C CALCULATE  THE  ELEVATION  FOR  EACH  CELL  IN  ONE  OUAORANT 


C and  The  same  value  can  he  A'^SIENFD  to  the  CORRAS- 
C SPCNDING  CELLS  IN  THE  OTHER  THREE  OUADRANTS 

c 

CO  ICC  I =1/  (NCELL+1) 

ROU  = (1-1)  • CELSIZ 
ZR?  = ROW**^ 

^ COL  = SQRT(ZRAO,?-;p?)  /FI.  OAT(CELSIZ)  + l. 

C 

C 

CO  ICO  CL=UNCnL 
COL  = CL-  1 

ZH  YP2=  ( COL  *C  ELS  I Z ) * •?  ♦ Z 
C 

C CALCLLATE  elevation  AT  CELL  USING  RIGHT  TRIANGiFS 

c 

ELEV  =SQRT(ZRA02  - ZHYP?) 

C 

c algorithm  to  place  value  in  corpessponding  cells 
C OF  all  four  quadrants 
c 

S = -1 

7C  S2=  -1 

71  C = S*COL  ♦ STCOL 

R = S2*(l-1)  ♦ STROW 

FL0(18/12/ARRAY(C»R))  = F L 0 ( ? A/ 1 ELE V ) 

FL 0 ( 30,6, ARR AY ( C ,R  ) ) = F L 0 ( 3 C , 6 , C A T ) 


1 F 

( S 

. £Q. 

1 ) 

GO 

TO 

pn 

S = 

1 

GO 

TO 

71 

8C 

IF  (S2 

.EC. 

S ) 

GO 

TO 

100 

T = S2 
S2  = S 
S = T 
GO  TO  71 
ire  CONTINUE 
RETURN 
END 

C SUBROUTINE  CYL<RAD,LFN) 

SUBROUTINE  CYL(RAD,LEN) 

IMPLICIT  INTEGER  (A-W) 

COMMON  BUF(?0),  A R R A Y ( S 00 , 7r ) , C F L S 1 Z , S T C OL , S T R OW , C A T , ANG 
C 

c 

PALFCEL  = CELSIZ/2 
ZRA02=  RAD**2 
NCOL  * LEN/CELSIZ 
NRCW  = RAD/CELSIZ 
C 
C 


CO  100  I 1= 1, (NROWA  1 ) 
1=  11-1 

ROW  = MCELSIZ 


i- 1 f >/  = s I3R  u / w A n , - row  • • M 
ru  1(0  j =1,NC0i. 

IJ  - j - 1 
C 

C 1- 1 A C F A (’  (>  R 0 F’  R I A T t ( L F V A 1 I o ' I ' f f 1 I '•  ' 0 I I ""i  ! '■  1 A ' ■*  v p 

C ANr  Ft  LOW  I HF  FF  NTF  RL  I M' 

C 


f L i)  ( 1 a , 

1 2 / AR  R A Y ( S r C 01 

+ J 1 » s 

1 R OW  < 1 

) ) 

= F 1. 

0 ( A L , 1 ? / r 1 F V ) 

T L C ( 1 8» 

1?,ARRAy(STroi 

♦ J ,1  , s 

r IJ  o w - ' 

' 1 ) 

= F 1 

ci(A/.  ,1  A/OlpV) 

F LC  ( m/ 

6/AKRAY  (SICT'L 

♦ J .1  » s 

I R-  ow  - 1 

[ ) ) 

s F 1 

D ( Tf  , ' , r A I ) 

FL  0 ( ?C/ 

6/AWRAYCSTCOI 

* J J 

I r 0 w ♦ 1 

’ > ) 

= F 1 

0 ( , C A T ) 

1 r t COM  I NUE 
R F TURN 
END 

C CL U ECU  TINE  IOUER(RAD/F|iV) 

SuEKCLTINE  TOWER(RAD/firV) 
itahlicit  integer  CA-W) 

rOM.VON  HUfCr’C)/  ARRAVC'  ' 0 , ? 0 ) , r f 1.  C 1 ' , S T C 01  , L T R OW  , C A T , A N G 
C 

C THE  SAKE  ALGCRIIU^  AS  LSED  fC  (.  HENISPhERF  IXCEpi 

C that  the  elevation  need  not  CALfilATAD  - IT  I'' 

C CONSTANT  value 
C 

t-At  f CEL  = CE  L S I Z /? 

NCELL  = FLCAI  (RAP)/FL OAicrriSIZ)  A.S 
AC2=  HA0«*^ 

ARRAyCSTCOl, STROM  = FIFV 
C 
C 

CO  1 (' L I = 1 / ( N C E L I.  + 1 ) 

ROW  = (1-1)  * CFLSIZ 
ZR2  = R0U**2 

NCOL  = SQRT(ZRAr2-Zr;>)  /FLn«T(fFtSI7)+i. 

C 
C 

CO  u:n  cL=i/NroL 

COL  = CL-1 

c 
c 

S = -1 
VC  S 2 = - 1 

71  C = S*CCL  ♦ STCOL 

R = S2*<1-1)  ♦ STROW 

f L D ( 1 e#  1 2/ AR  R A Y ( C/ R ) ) = f L f1  ( 2 A # 1 ' » f L > V ) 
FLU(30/6/ARRAY(C,R))  = 

If  ( S . E Q . 1 ) GO  TO  R r 
S =1 

ro  TO  71 

C If(S2  .EQ.  G)  GC  TO  mo 


C? 


( 


ivrw’^N'l 


CO  TO  n 
t C 0 N T 1 N U t 
RETURN 
END 

SUbROUTlNE  CHKRPO 
SUbROUTINE  CHKRbC 
IMPLICIT  INTRf3FR  (A-W) 

COMMON  BUE(20)#ARRAY(5nn,7r)/CEL»:!Z/STCCL/!?TR0W 

CO  icr  1*1/20/4 

CATEG  = I/4+1 
CO  100  J=1/4 
CO  ICO  <=1/4 
fiOU  =J  ♦ STROW 
COL  = l^K+STCOL 
R0Vj2  = I+K+STROW 
C0L2  = J+STCOL 
ARRAY{COL/ROW)  = CATFr, 

C ARRAY(C0L2/R0W?)=CATEG 
RETURN 
END 

subroutine  FILLIN(LC/BC/LR/RR/CATFG) 

SUBRCLTINE  FILLINTLC/BC/LR/RR/CATFG) 

THIS  SUBROUTINE  CHECKS  ROTATED  SURFACES  FOR  'HOLES 
AND  FILLS  THEM  IN  WITH  APPROPRIATF  FLFVATION  AND 
CATEGORY  IF  IT  FINOS  ANY 

implicit  INTEGFR(A-W) 

COMMON  BUF (20) / ARRAY(SnO/7r) 

DATA  HUN/0100/ 


BEGIN  AT  LEAST  ROW  AND  COLUMN  AND  PROCFFO  TO  BIGGFST 
ROW  AND  COLUMN  TO  CHECK  FOR  HOLES 


DO  ICO  I»LR/BR 
CO  ICC  J * L C/BC 

IF  TEE  CELL  IS  N0N7FP0  IT  HAS  HFEN  ASSIGNED  A VALUE 
AND  IS  OK.  PROCEED 

IF  ( ARRAY< JE 1 / I ) .NE.  0)  GO  TO  100 

THIS  CELL  IS  ZERO.  IF  BOTH  OF  ITS  NFIGHOORS  (TO  LEFT  AND 
RIGHT)  ARE  NON  ZERO/  THEN  WE  HAVF  FOIlfgD  A HOLF 
otherwise  THIS  CELL  IS  NOT  FART  CF  THE  OBJECT 
AND  can  BE  IGNORED 


I F ( ARRA  V ( J /I  ) .FQ.r  .OR.  ARRAY(J*?,I)  .FT.  n)  r,0  TO  100 


CRINC  LABORATORIES 


Chemical  Engineering  Low  Temperature  Laboratory 

Remote  Sensing  Laboratory 

Flight  Research  Laboratory 

Chemical  Engineering  Heat  Transfer  Laboratory 

Nuclear  Engineering  Laboratory 

Environmental  Health  Engineering  Laboratory 

Information  Processing  Laboratory 

Water  Resources  Institute 

Technical  Transfer  Laboratory 

Air  Pollution  Laboratory 

Satellite  Applications  Laboratory 


